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ABSTRACT 



When a beam of extraordinarily polarized light is directed through a poled, single 
crystal of BaTiO^ it fans in the direction of the optic axis. This thesis research 
investigates the origin of this fanning, or asymmetric self-defocusing, in barium 
titanate. Experimental observations and theoretical simulation suggest that this 
phenomenon is due to near forward stimulated photorefractive scattering. 
Additionally, observations of unique variegated beam patterns during phase 
conjugation are reported, including internally pumped rings and motion invariant 
patterns. 
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I. INTRODUCTION 



During the last decade the photorefractive crystal barium titanate (BaTiO^) has 
received a large amount of attention within the scientific community. Barium titanate 
has long been of interest to acousticians and solid state physicists because of it's 
photoacoustic and piezoelectric properties, but only recently have optical physicists 
become very interested in this crystal. The cause of this sudden and intense interest 
was the discovery of photorefraction and optical phase conjugation. Originally found 
in liquids and gases [Ref 1], optical phase conjugation has the ability to reconstruct a 
distorted wavefront (sometimes loosely referred to as time reversal). Barium titanate is 
one of many photorefractive crystals which can exhibit phase conjugation. 

There exist a number of intriguing optical properties of single crystals of barium 
titanate, but by far the most interesting is its ability to phase-conjugate at continuous 
wave low powers, with no external apparatus. This effect, known as self-pumped phase 
conjugation (SPPC), is the subject of extensive current investigation. The goal of this 
research is to determine the origin of self pumped phase conjugation in BaTiO^. 

Currently there exist two explanations of the cause of SPPC in BaTiO^. One 
theory, developed by Feinberg and his students at UCLA [Ref. 2], attributes the 
production of the phase conjugate reflection to the photorefractive effect via 
Degenerative Four Wave Mixing (DFWM). According to this theory, an asymmetric 
change in the index of refraction across the path of an incident beam causes the beam 
to bend. This bending eventually leads to a geometry which, through two internal 
reflections (often termed a CAT corner), allows the beam to intersect itself, leading to 
DFWM. The second theoiy^ originating with Lam and his co-workers at Hughes 
Research Laboratory [Ref 3], attributes SPPC to Simulated Photorefractive Scattering 
(SPS), or Two Wave Mixing, similar to that exhibited in gases and liquids with two 
frequency disparate input beams of very high powers. With the goal of proMng one or 
the other of these theories, this thesis begins with a general background and theory of 
phase conjugation as it applies to BaTiO^ (Chapter II). Chapter III describes PSI.M 
(Photorefractive SIMulation), a simulation of the phenomenon of beam fanning due to 
a change in index of refraction within the incident beam path arising from the 
photorefractive effect. This simulation enables the prediction of the change in index of 
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refraction and the index profile necessary to account for SPPC with only the 
photorefractive effect taken into account. 

Chapter IV introduces experiments to determine the effect of incident beam size 
on various crystals of BaTiO^. Chapter V summarizes the results and analyzes the 
data presented. The conclusion is reached that Stimulated Photorefractive Scattering is 
responsible for the process of beam fanning which leads to self pumped phase 
conjugation in BaTiO^. The actual origin of the phase conjugate beam is attributed to 
both SPS and DFWM. 

Additionally, Appendix A contains a summary of the work I did at Los Alamos 
National Laboratory under the direction of Dr. R. A. Fisher in February and March, 
1986. This work describes and catalogs internal beam production in barium titanate 
during SPPC. Some internal beam patterns not previously reported are introduced and 
explanations for these phenomena are proposed. The results of the experimental work 
in the body of this thesis have been submitted for publication to Optics Letters. An 
expanded version of Appendix A was presented at the XIV International Quantum 
Electronics Conference {1986), and has been submitted for publication in the Journal of 
the Optical Society of America B in collaboration with R. A Fisher and A. V. Nowak. 
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II. BACKGROUND AND THEORY 



A. INTRODUCTION 

Optical phase conjugation, first discovered in 1972 [Ref. 1], relies on a specific 
nonlinear property of a substance to create the formation of the complex conjugate of 
an incident wave. The reflection involved is not the familiar specular reflection seen in 
everyday life, but a reflection that retraces the path of the incident beam exactly, 
regardless of the direction or angle of incidence. The exact retracing of the incident 
beam path and the formation of the complex conjugate of the incident wavefront are 
the two unique and most important effects that make phase conjugation both 
interesting and useful. In barium titanate this phase conjugate beam is directly a result 
of the photorefractive effect. 

B. ORIGINS OF THE PHASE-CONJUGATE BEAM 

The origin of the phase conjugate beam in photorefractors is light which is Bragg 
scattered by an index grating formed by the light in the material. This process is 
known as wave mixing. There exists both two-wave mixing (scattering) and four-wave 
mixing in barium titanate. Both of these processes will be treated theoretically in 
general terms; the specifics of wave mixing in BaTi 03 will be treated later. 

To understand the origins of phase conjugation one must begin with the wave 
equation for a homogeneous, isotropic, time-independent, nonmagnetic and 
nonconducting dielectric. The polarization term is divided into its linear (Pj) and 
nonlinear (P^j) parts, and the linear portion is absorbed into the linear permittivity of 
the medium, £. This leaves the wave equation as: 

, d^E d^P. 

v2e-cu^= (2.1) 

Ol~ Ol~ 

For ease of manipulation an approximation to this equation is made that reduces 
the second order nonlinear wave equation to a first order equation.^ This 
approximation is derived by Fisher and Yariv [Ref 5; pp. 9-11] and for it to be vahd, 
the envelope of the pulse must not change appreciably during an optical period, hence 

^Some calculations can be made without this approximation but they are quite 
tedious (e.g. [Ref 4]). 
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the term Slowly Varying Envelope Approximation. In making the approximation it is 
necessary to assume that the electromagnetic waves under consideration are plane 
waves, and also to extract only the portion of the nonlinear polarization that 
propagates with the same frequency and wave vector as the modulation function of the 
electric field. This, referred to as the phase matched portion of the polarization, is 
given by 



P(z,t) exp±i(cot-kz), 
in the wave given by 



( 2 . 2 ) 



E = E(z,t) exp±i(o)t — kz). 



(2.3) 



Here the boldface letters (E and P) are the fields and the conventional letters (E and P) 
are the field envelopes, which vary little in an optical period and in an optical 
wavelength. All other terms are ignored because they are not phase matched and 
therefore cannot couple to the electric field. The Slowly Varying Envelope 
Approximation, or SVEA, is then invoked. That is: 



, dE d^E 

Ik^EI >>lk— I>>!-3- 
dz oz^ 



(2.4) 



Fisher and Yariv then show that the wave equation with SVEA becomes a first-order 
equation relating the envelope functions E and P: 



dE , dE i(oVl> 
— + = ±- -r- P. 

dz ^ dt 2 



(2.5) 



This result is the SVE.A wave equation for a plane wave traveling in the — z direction, 
which shows that, given a nonlinear polarization, an electric field can be produced. 
This development can be generalized to the case of many waves present within a single 
medium, resulting in the origin of the nonlinear processes that can lead to phase 
conjugation. 

The polarization term can often be written as a power series expansion: 



12 



P = EX(E) = X^^)e + x(2)e^ + x(^)e^ + • • • 



( 2 . 6 ) 



where X(E) is the susceptibility of the material in which the wave is traveling. The 
expansion of the susceptibility in a power series separates the effects of the nonlinear 
material into distinct groups that categorize nonlinear effects [Ref 5; p. 13]: 

• X^^^ - These are the linear properties that are the subject of classical optics and 
are not directly responsible for nonlinear effects. 

• X^^^ - The second order effects are in general called three-wave mixing and 
include: second-harmonic generation, optical rectification, parametric mixing 
and the Pockels effect. 

• X^^^ - Third order effects include some of the most popular effects responsible 
for phase conjugation such as: third-harmonic generation, nondegenerate four- 
wave mixing, Raman scattering, ac Kerr effect (degenerate four- wave mixing), 
dc Kerr effect, Brillouin scattering, two photon absorption, dc-induced 
harmonic generation and Stimulated Photorefractive Scattering. 

The two methods of creating a phase-conjugate reflection that are directly applicable to 

the experimentation described here are scattering and degenerate four wave mixing. 

1. Four-Wave Mixing 

In 1977 Hellwarth [Ref 6] showed that it was possible to generate a time- 
reversed repUca of any monochromatic-beam wave pattern. The method employed the 
interaction of an incident beam with counter-propagating pump waves in a 
homogeneous, transparent, nonlinear medium (see Figure 2.1). This method, called 
four-wave mixing, has become one of the most popular methods for phase conjugation 
because of the almost universal availability of the necessary equipment and proper 
media. 

The description of wave mixing begins with the nonlinear polarization derived 
from the energy function. As shown by Pepper and Yariv [Ref 5: pp. 26-33], isolating 
only the third order term gives 

- «'l2pcElE:Ep ■ (2-7) 

I consider here degenerate four-wave mixing (DFWM), where 

(0j = C 02 = (0p = (0^ represents frequencies of the four beams. The subscripts 1 and 2 
refer to two strong, precisely counterpropagating pump beams, the subscript p refers to 
the weak probe beam, and the subscript c refers to the phase-conjugate output wave. 
It is assumed that the following conditions are met: 
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k. 



‘2 



Figure 2.1 Schematic of the process of four-wave mixing [Ref 5; p.421], 

a) intense (nondcpletcd) counterpropagating planewave pump beams of equal 
intensity (kj, k-,), 

b) a weak probe beam (kp), 

c) both the intense pump beams and the weak probe beam are at the same 
frequency (O. 

Since the two pump waves are very strong compared to the probe beam they form a 
standing wave. To obtain the steady-state solutions, we can set 



So. assuming the SVEA and considering only the phase matched portion of tlie phase 
conjugate wave, then from Equation 2.5 




( 2 . 8 ) 




(2.9) 
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The nondestructive buildup of the electric field amplitude can occur only if the product 
of the exponential phase factors in Equation 2.9 and is small (on the order of zero). 
If this is not so, then all of the successive contributions will not be in phase to add 
constructively. This total phase factor is: 

exp[-i(kj + l4-kp.k^)»r]. (2.10) 

Since waves 1 and 2 are counterpropagating then, 

kj + kj = 0. (2.11) 

Although in practice this may not be exactly so because of nonlinear phase shifts due 
to unequal pump-wave intensities or non-aligned beams, it is very close [Ref. 5: p. 30]. 
Therefore the phase factor can be zero only if 




( 2 . 12 ) 



Here it is evident that the nonlinearly generated field propagates in the opposite 
direction of the probe wave. 

Assuming that the factor in the exponential of Equation 2.10 is equal to zero, 
then combining Equation 2.7 and Equation 2.9 yields (within a multiplicative constant) 
[Ref. 5; p.36] 



dE 

dz 



■^= ico 






(2.13) 



This equation was solved concurrently by Yariv and Pepper [Ref. 7] and 
Bloom and Bjorklund [Ref. 8], and the solution was shown to be: 



Ep(z) = 
Ec(z) = 



. ik| sin(iklz) cos[lkl(z-L)| 

-1 E (L) + E (0) , 

ik| cos(lkiL) cos(ik|L) P 

cos(|k|z) ^ ^ * sin[|k|(z-L)] ^ 

E.,(L) + ik E„ (0). 

cos(ik|L) |k|cos(|k|L) P 



(2.14) 

(2.15) 
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Where 



k* = (0 X^EjE2 , 



Ve 



(2.16) 



E^(L) is the amplitude of the conjugate wave at its origin (the phase conjugate mirror), 
and Ep(0) is the amplitude of the probe wave at its origin. 

In practice, the amplitude of the phase conjugate wave at its origin is zero (i.e. 
E^(L) = 0), so that at the input point (z = 0) of the probe beam 



E^(0) = i — tan(|k|l)E >). 
|K| ^ 



(2.17) 



Equation 2.17 shows that at z = 0 (the origin of the probe beam) the reflected field is 
proportional to the complex conjugate of the incident field. Therefore there exists a 
phase conjugate reflection. 

2. Scattering 

Phase conjugation by stimulated scattering was first reported by Zel'dovich, 
Popovichev, Ragul'ski and Faizullov in 1972 [Ref 1]. As stated by Hellwarth 
(Ref 5: p. 171]: 

In essence, the nonlinear polarization density, mediated by some driven (Raman-, 
Brillouin-, etc.) active excitation of the optical medium, couples normal-mode 
solutions of the linear Maxwell equation so as to create a new set of normal 
modes in backscattering with, in effect, complex propagation constants 
representing index change and gain. 

For example Stimulated Brillouin Scattering, or SBS, is a direct result of the 
electrostrictive effect seen in some materials. This effect is seen when the polarizability 
of a material is a function of pressure. In these materials it is possible for light to be 
scattered by a sound wave allowing che coupling of a pair of light waves to a pressure 
wave. Another type of scattering sometimes responsible for phase conjugation is 
Stimulated Raman Scattering or SRS. SRS is the scattering of light from molecular 
vibrations or rotations in a gas or liquid, or longitudinal-optical phonons in a solid. In 
SRS light waves are coupled because of polarizability changes with a molecular 
coordinate. Of particular importance to this research is Stimulated Photorefractive 
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Scattering in barium titanate. In SPS it is the photorefractive effect which is 
responsible for the production of a periodic index of refraction grating; this grating 
then Bragg scatters light in the backward direction which comprises the phase 
conjugate beam. 

The mathematics involved in all three types of scattering is similar and can be 
treated together . To this end I outline the work of Hellwarth [Ref 5: pp. 177-180] and 
discuss the existence of the phase conjugate wave assuming scalar or symmetric 
scattering with a monochromatic, multimode, incident wave in a wave guide. The 
results are largely independent of the nature of the wave guide provided that a specific 
TEM mode that propagates in the forward direction also propagates in the backward 
direction. 

A wave E^(r) may be expressed as: 



where P is the attenuation coefiicient, is the normalized transverse-mode pattern, 
and the propagation constants may vary with wave angular frequency, that is k^ = 



The electric field of the backscattered wave may then be written as 
[Ref 5: p.l80] 



where k^^ = kj^(co), and y is used to balance any nonlinear terms from the fields 
interaction with the medium. If y = 0 in the above equation, then this is the solution 
to the linear Maxwell equations for a wave in a wave guide. 

If in Equation 2.19 y5=0. then there is one configuration that allows the y of 
one of these modes to have a positive reai part significantly greater than that of any 
other. This means that a single wave (refiection) may dominate the backscattered 
waves if the necessary conditions are met. The desired result is that 




(2.18) 




^2 "" ^ - yz/2] , 



(2.19) 



E^{t) = nEi'‘(r), 



( 2 . 20 ) 
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where n is a constant. Substituting Equations 2.18 and 2.19 into Equation 2.20 gives 
the polarization density necessary to balance the extra terms in Maxwell's equations 
that come from the last term in Equation 2.19 . According to Hellwarth (Ref. 9] the 
necessary conditions are satisfied if 






4 f 



m 



*4. *6 K^:- A.A.B 
ij n niijn 1 j n 



( 2 . 21 ) 



where 



^mijn = exp[AkLz]dz, (2.22) 

and 

Ak = kJco) + ki(v)-kj.(v)-k„((o). (2.23) 

In Equation 2.21 I have assumed that the scattering considered is the pure scalar or 
symmetric type, reducing the fourth rank susceptibility tensor to a scalar quantity. 
This quantity, and all other constants, I have absorbed into i]. 

The necessary conditions may be determined for special cases, either by exact 
solution or perturbation methods developed by Hellwarth. These solutions indicate 
that a backscattered solution that is nearly phase conjugate predominates because it 
has gain nearly twice that of any other mode, provided that [Ref 5: p. 177]: 

a. the interaction length is not too long, 

b. the Stokes shift (v — to) is not too large, 

c. the number of guide modes falls within a certain large range, and 

d. the total number of guided modes is not too large. 

The solutions that led to these conditions, presented by Hellwarth in 1978 [Ref 9], 
unfortunately do not take into account competing nonlinear e.Tects, and thereibre are 
good only to the extent that they predict the major effects observ'ed in experimentation. 
It is observed that the backscattered wave is not always an exact phase conjugate of 
the input wave, possibly due to some of these competing effects. These analyses apply 
equally as well to SPS, SBS and SRS [Ref 3j. 
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C. DISTORTION CORRECTION BY PHASE CONJUGATION 

Although not immediately applicable to this investigation, the ability of a phase 
conjugator to correct for phase distortions is one of the primary aspects driving 
research in this area. As a practical point, this ability is an easy way to determine that 
any given observed beam is a phase conjugate of another. 

It was shown above that, at least for the two cases considered, a wave can be 
generated propagating directionally opposite to an input wave, with it's envelope 
function proportional to the complex conjugate of the input beam's envelope function. 
This phase conjugate wave may be considered a separate wave that must satisfy the 
same wave equation as the input wave. This means that if the wavefront of the input 
wave is distorted in some way, then the phase conjugate wave will propagate 
backwards along the same path and its envelope function will remain everywhere the 
complex conjugate of the input wave's; that is, its wave fronts will coincide with the 
input wave fronts at every point. Thus a distorted wave front will be restored to its 
original shape as the phase conjugate wave traverses the distorting medium in the 
reverse direction. The proof is quite simple, as shown by Fisher and Yariv [Ref 5: p. 
17], and is accomplished by complex conjugating the wave equation for the forward 
going wave. The result is the equation for the conjugate wave propagating in the 
opposite direction. So a wave traveling in the opposite direction of the incident wave, 
with an envelope function the complex conjugate of such a wave, satisfies the same 
wave equation as the incident wave. This means that the backward going wave will 
remain everywhere the complex conjugate of the incident wave, and therefore a 
distorted plane wave will be restored to its original configuration when it is passed back 
along the same distorting path. This distortion correction capability is the most visible 
and fascinating effect attributable to phase conjugation at the present time. 

An incomplete list of some applications of phase conjugation with emphasis on 
those currently using barium titanate as the phase-conjugator is [Ref 5,10,11]: 

a. Holography 

b. Distortion correction 

c. Laser oscillators 

d. Interferometry 

e. Laser beamsteering 

f Associative memory investigations 

g. Image amplification 

h. Laser coupling 
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i. Image addition, subtraction and enhancement 

j. Optical bistability 

k. Automatic pointing and tracking. 

D. PHOTOREFRACTORS AND FOUR WAVE MIXING 

In photorefractive materials the index of refraction changes with incident optical 
energy. This effect was first discovered in the 1960's [Ref 12,13] and has since led to 
the development of optical recording sensitivities comparable to that of the silver halide 
emulsions used in photography [Ref 14]. The process involved in the effect is 
summarized by Feinberg [Ref 5: p. 418] as follows: 

a. Light causes charge to migrate and separate in a crystalline material. 

b. The separation of charge produces a strong electrostatic field. 

c. The electrostatic field causes a change in the refractive index of the crystal by 
the linear electro-optic effect.^ 

In the photorefractive effect the important factor in the change in index of 
refraction is not the intensity of the incident light, but the total energy of the incident 
radiation; the power only determines the speed of the reaction. Also because of the 
type of effect concerned, the actual result is dependent on the relative intensity of the 
incident radiation. 

Consider a thin beam of light incident on a photorefractive crystal that is 
otherwise dark. The light will liberate charge carriers from the lighted area of the 
crystal and these carriers will then settle back into a place within the crystal structure. 
If the carriers settle back in the area of the crystal that is lighted they will again be 
liberated. If, however, the carriers land in a dark area of the crystal, they will stay 
where they land due to the small dark conductivity of the crystal. Eventually there will 
be few charge carriers in the lighted area and many in the dark area. This 
disproportionate distribution creates an electric field between the dark area and the 
lighted area that creates a change in the index of refraction of the crystal through the 
linear electro-optic effect. 

This process of migration has been aptly described by R. A. Fisher as the 
ilyswatter effect [Ref. 16]. It is as if the charge carriers were flies at a picnic table and 
the light was a flyswatter suspended from above, slowly swinging across the table top. 

^The linear electro-optic effect is the change in the index of refraction caused by 
(and proportional to) an applied dc electric field. This effect is seen in crystals that do 
not exhibit inversion symmetry, such as barium titanate in the tetragonal phase 
[Ref 15: p. 275]. 
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Each time a fly lands on the table within the area of the flyswatter's reach, he will 
immediately fly up as the swatter returns toward him. If, however, the fly lands on the 
table in an area away from the swatter he will be left alone and will stay on the potato 
salad. In time therefore, there will be no flies on the table in the area where the 
flyswatter can reach and many where it cannot. This process effectively pumps the 
flies to regions on the table which are away from flyswatter. 

The model of charge transport described above is known as the hopping model 
and was first proposed by Chen to explain the change in refractive indices in LiNbO^ 
and LiTaOj [Ref 17]. There also exists a diffusion model proposed by Amodei 
[Ref. 18] that may dominate when the spacing between light and dark areas within the 
crystal is less than one micrometer. 

From the flyswatter analogy it is clear that what is important is the optical 
intensity relative to the ambient. This entire explanation assumes the existence of 
charge carriers in a crystalline material. In most materials exhibiting the 
photorefractive effect the origin of the charges is unknown, but it is assumed that they 
inhabit low-lying traps formed by impurity sites in the crystal. For example, in 
BaTiOj, Fe^"*" and Fe^**" impurities may create the charge carriers as shown by Klein 
and Schwartz [Ref 19]. It has also been proposed, and there is extensive evidence to 
support the theory, that oxygen vacancies are the donors and barium vacancies are the 
acceptors [Ref 20]. Whatever the dominant species is, it is assumed that there are a 
large number of empty receptor sites in the material available to the liberated charge 
carriers. 

If instead of a single beam, two beams of identical wavelength and equal 
intensities (one called a pump wave and one called a probe wave) intersect within a 
photorefractive crystal, the interference pattern they form will cause a periodic optical 
field to be set up in the region of intersection, resulting in a periodic charge separation. 
Through the electro-optic effect, a periodic variation in index of refraction will result. 
This index grating is then able to Bragg scatter a pump wave in the opposite direction 
to produce a phase conjugate beam. The formation of the index grating is shown in 
Figure 2.2, adapted from Feinberg [Ref 5: p. -122]. 

It is important to note that photorefractive phenomena are not merely four- wave 
mixing phenomena. Four- wave mixing is a consequence of an index change at some 
given point at some given time (x,t). Photorefractive phenomena however, are 
dependent on the history of the medium and are not spatially local effects. That is, the 
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Figure 2.2 Index grating formation due to the photorefractive cftect [Ref. 5). 

effects at (x,t) are a consequence of actions at (x- 6 x, t- 6 t). This is clearly seen in 
Figure 2.2. 

The top curve of Figure 2.2 shows the intensity as a function of the position 
along the path of the beam. The second curve is the charge density due to this 
intensity variation. The electric field intensity is presented in the third curve (note the 
important V 4 pei'iod shift relative to the charge). Finally the variation of the actual 
index of refraction within the cr>'stal is shown in the last curve. 

The spatial shift between the change in index of refraction and the variation in 
intensity due to the interference pattern (one fourth of the grating period) is an 
important effect. This happens because the peak of the charge distribution is a center 
of right-left symmetr}', and therefore there can be no electric field at that point. This 
shift leads to the ability to transfer power from one beam to another, a process known 
as two-beam coupling. Two-beam coupling, first reported by Staebler and Arilodei 
[Ref 21], is due to the scattering of one beam into another by Bragg diffraction as they 
pass through the grating they have mutually produced. The result is that one beam 
experiences gain at the expense of the other. 
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One interesting efTect of four-wave mixing in photorefractors is seen in the phase 
conjugate mirror with continuous wave gain. This was first demonstrated by Feinberg 
in 1980 [Ref. 22] using the configuration shown in Figure 2.3 . 



Figure 2.3 Formation of an optical cavity with a phase conjugate mirror [Ref. 5). 

In this configuration a phase conjugate mirror (PCM) and an ordinaiy mirror (Ml) 
can form an optical cavity in which a beam will oscillate. One of the more dramatic 
demonstrations of this type of oscillation is reported by Feinberg [Ref. 5: p. 436] in a 
configuration using an ordinary kitchen spatula as one end of the cavity. Just as in 
any conventional laser oscillator, the oscillations can begin with noise (as little as a 
single photon) that makes one successful round trip within the cavity. With each 
additional round trip the system sees gain due to two-boam coupling and in a short 
period of time a visible CW beam may be observed between the mirror and the phasc- 
conjugator. 

E. STIMULATED PHOTOREFRACTIVE SCATTERING 

The second theory which accounts for phase conjugation in BaTiO^ was 
originally proposed by Lam in 1985 [Ref. 3|. Lam proposed that the origin of tiie 
phase conjugate beam is two wave mixing or scattering. This process of SPS begins 
with random inhomogeneities in the photorefractive medium which give rise to noise 
photons from the Rayleigh scattering of the input beam. In the words of Lam [Ref 3j: 

The noise photons propagating along with [sic] axis of the input beam undergo 
parametric scattering, resulting in an unidirectional gain. The phase conjugate 
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wave arises from the component of the scattered wave which experiences the 
maximum gain. 

Later the same year Chang and Hellwarth conclusively demonstrated SPS in BaTiO^ 
[Ref 23]. To do this they had to index match the surfaces of the crystal to make 
accessable angles which would not produce the characteristic corner reflection (i.e. the 
CAT corner) producing DFWM. 

Later, Valley [Ref 24] deflned SPS as the selective amplification through the 
photorefractive effect of optical radiation scattered by medium inhomogeneities. Thus, 
SPS is a two-wave mixing process in which the incident beam interferes with scattered 
light to create the necessary grating through the photorefractive effect. As soon as the 
gain for an individual noise beam exceeds its loss, then that beam will become 
dominant, increasing the grating strength and hence the backscattered beam gain. This 
process will continue until a steady state value is reached. 

The mathematical equations, as derived by Lam, are of the identical form to 
those encountered in stimulated scattering in a waveguide (presented above), so the 
analysis of phase conjugation by SPS has already been addressed. 

F. BARIUM TITANATE 

1. The Crystal 

The experiments reported here were performed with one of the more popular 
photorefractive crystals capable of creating a phase conjugate reflection. Above the 
Curie temperature of 120°C, BaTiOj is a cubic crystal and therefore centrosymmetrical 
and non-photorefractive. Between 5°C and 120®C however, the crystal is in the polar 
phase and has a tetragonal symmetry. The axis of fourfold rotation is the optical axis 
or c-axis. In the tetragonal phase barium titanate is ferroelectric, photorefractive, 
photoconducting, birefringent, electro-optic, pyroelectric, piezoelectric, photoacoustic 
and (slightly) electro strictive. The many possible effects make the crystal itself very 
difficult to understand, since often the results of the various effects cannot be 
separated. All results reported here involve BaTiO, in the tetragonal phase. 

Since BaTiO, is cheinicallv and mechanicallv stable at room temoerature. it is 
used extensively in many capacities and is therefore well studied. Jona and Shirane 
[Ref 25: pp. 108-215] report extensively on the structure, properties and the effects of 
temperature, stress and electric fields on barium titanate in the tetragonal phase. 
Although at the time of printing the photorefractive capabilities were unknown, their 
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data represent a comprehensive picture of the crystal and forms the basis for analysis 
of the solid state aspects of barium titanate. More recent but less comprehensive data, 
emphasizing the optical aspects of BaTiO^, have been compiled by Wemple, et. al. 
[Ref 26]. 

When used as a photorefractor, it is necessary that a single crystal of BaTiO^. 
be used and that the crystal be poled into a single domain. Poling is accomplished by 
cooling the crystal through the Curie point while an electric field on the order of 10^ 
V/cm is applied across the optical axis. This poling of the crystal results in an electro- 
optic tensor which has a very large [42] element. Since the photorefractive effect uses 
the linear electro-optic effect, it is important to maximize this [42] element. This 
maximization is accomplished by using extraordinarily polarized light incident at an 
angle between zero and 90° to the c-axis. 

2. Beam Fanning 

When a beam of extraordinarily polarized light is incident on a poled, single 
crystal of barium titanate, it is seen to asymmetrically defocus in the direction of the c- 
axis prior to the onset of phase conjugation. This process, called beam fanning, was 
first discovered by Feinberg in 1981 [Ref 27]. Beam fanning is clearly visible in Figure 
2.4 (the optic-axis is designated by the arrow). This fanning can be explained by either 
an asymmetric change in index of refraction within the incident beam path, or by 
stimulated scattering. 

a. Deflection by Photorefraction 

Beam fanning was proposed to be due to an asymmetric change in the 
index of refraction caused by the photorefractive effect. This process is represented in 
the curves of Figure 2.5 (from Feinberg [Ref 27] ). The top curve shows the intensity 
curve of a Gaussian beam. The second curve shows the net charge due to 
displacement of charge carriers by the incident light. The third curve shows the 
induced electric field, and the last curve shows the asymmetric change in the index of 
refraction. It is assumed that at the center, and most intense portion of the beam, the 
gradient of the index of refraction is greatest. Therefore the center portion of the beam 
creates the observed effect, and the fanning is induced beginning at that point. In 
other words, beam fanning begins in the center of the beam. This is one of the major 
discriminating factors and will be referred to later.^ The second major point to consider 

^Note that even if the intensity profile of the incident beam is not smooth, the 
top curve of Figure 2.5 represents the envelope that the beam intensity must fall 
within. 
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Figure 2.4 BaTiO^ exhibiting SPPC. 

about this theory is that as the diameter of the beam decreases the electric field 
increases, since the linear distance between charge distributions decreases. Since the 
photorefractive effect is proportional to the electric field, then beam fanning must 
increase with a decrease in beam diameter. This is clearly seen by examining the 
equations describing the electric field and the subsequent refractive index change. 

The electrostatic field induced by an incident beam in a photorefractor is 
given by [Ref. 27]: 



-2AdT(xx + vy) T T 1 

E(x) = exp[(x2 + y2)/o) 2], (2.24) 

where e is 'he charge of the charge carrier, is Boltzmann's constant, T is the 
temperature of the ciystal, and (0^ is the incident beam diameter. The change in index 
of refraction in a crystal of BaTiO^ for an extraordinary ray due to this field is 
approximately given by 
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Figure 2.5 The process of beam farming due to photorefraction 
for a single incident beam [Ref 27J. 

n(y) = -|E(x)*y] -^-^^(2n^^ng-r^2sin20), (2.25) 

where n(0) = (n^" cos0 + sin0)’^^^^ , 0 is the angle the incident beam makes with 

the non c-axis of the ciystal, n^ is the index of refraction for an ordinary ray, n^^ is the 
index of refraction of an extraordinary ray, and all terms that do not include the r^, 
element of the electro-optic tensor have been dropped, since its value exceeds all other 
elements by an order of magnitude. 

Inspection of these equations clearly shows that the change in index of 
refraction is inversely proportional to the diameter of the incident beam. 
b, Near Forward-SPS 

In 1985 Lam proposed that beam fanning is due to an effective nonlinear 
index of refraction which causes the cry'stal to behave like a diverging lens [Ref 3]. 
This nonlinear index of refraction (~10'^) is due to an intensity dependent frequency 
shift of the incident beam due to the bulk photovoltaic effect. Recently Valley 
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(Ref. 24] has suggested that beam fanning is actually near forward scattering due to 
SPS, in an extension of Lam's theory. 

As in other types of scattering, SPS is the selective amplification of 
radiation which is scattered by inhomogeneities in the medium. The amplification is 
accomplished through the photorefractive effect. This is schematically diagrammed in 
Figure 2.6. 




Figure 2.6 The process of SPS [Ref 24]. 

In a crystal where one charge carrier dominates, the gain-length product 
per absorbed energy per volume at turn on is given by [Ref 24] 



G 



to 



2n^[£Cjl + E,/EJ] ’ 



(2.26) 



where r^j^ is the eifective Pockeis coerncient for extraordinarily polarized light [Ref 2Sj, 
k = 2rt/X, e is the carrier charge, T is the space charge field and refractive index grating 
decay rate, £ is the dielectric constant (which is anisotropic), and 
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where Ag is Boltzmann's constant, ^ is the carrier mobility (which is also anisotropic) 
and T is the temperature of the crystal. Finally, Ig^ is the effective interaction length 
which Valley models with 

Igfy = L{1 - exp[-<o^/(L sinG)]}, . (2.2S) 

where L is the length of the crystal in the direction bisecting the scattered and incident 
beam, co^ is the incident beam diameter, and 

e = (e, - 8p)/2 . (2.29) 

The important item to note is that the gain-length product is proportional 
to the effective interaction length given by Equation 2.28 . Therefore the amount of 
scattering due to SPS (both forward and backward) is exponentially dependent on the 
incident beam diameter in such a way that an increase in beam diameter will increase 
the scattering, and a decrease in beam diameter will decrease the scattering; this is 
opposite from that expected by the theory discussed in Section F.2.a above. Also, 
since SPS is dependent on the interaction of random noise photons with the incident 
beam, one would expect that beam fanning (near forward-scattering) could commence 
at the periphery of the beam in addition to the center. This scattering of energy should 
be more obvious at the edges of the incident beam since the intensity is lower there. 
Also, any light scattered toward the c-axis from inside the beam would take power 
from the periphery through two beam coupling, thus adding to the loss of intensity 
toward the beam edge. Considering these points, if the intensity profile of the incident 
beam were to be monitored during the fanning process, one should expect to see the 
intensity diminish from the periphery toward the center. Again this is opposite the 
effect expected from the prior theory. 

3. Self Pumped Phase Conjugation 

Beam fanning is intricately involved with seif-pumped phase conjugation 
(SPPC). SPPC in BafiOj was iirst discovered by Feinberg soon after the discovery of 
beam fanning [Ref. 2]. According to Feinberg, the process of beam fanning (due to an 
asymmetric change in index of refraction within the beam path) may result in a portion 
of the incident beam being bent in such a way that it is internally reflected off the 
corner of the crystal (see Figure 2.4). This reflected beam then intersects the incident 
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beam causing a steady state interference pattern. The light and dark areas of the 
interference pattern cause an index grating (through the photorefractive effect) which 
then Bragg scatters the incident beam in the process of DFWM, creating a phase 
conjugate beam. Lam alternatively theorized that the origin of the phase conjugate 
beam is two-wave mixing (SPS) as described above, and that auxiliary' beams were not 
essential to the process. 

Of interest also is the time for phase conjugation to begin in barium titanate. 
The time to produce self-pumped phase conjugation can vary from a few seconds to 
tens of minutes depending on the intensity of the incident beam and the configuration 
used. Because of the typically large time for the onset of phase conjugation in BaTiO^, 
there’ is extensive work presently underway to try to decrease it. Some work on 
changing crystals after growth has been reported [Ref 20], but the present emphasis is 
on the doping of the crystals during growth [Ref 29]. 

The exact processes involved in beam fanning and phase conjugation are still 
unresolved in the literature. In the following chapters I will examine these theories as 
they apply to observed effects and, since in the case of beam fanning there are 
mutually exclusive indications, I will show that beam fanning must be due solely to 
near forward SPS. 
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III. PSIM: A SIMULATION OF THE PHOTOREFRACTIVE EFFECT 

IN BARIUM TITANATE 



A. INTRODUCTION 

To investigate the photorefractive effect in barium titanate, I designed a 
computer simulation which plotted the internal beam pattern within the crystal. This 
simulation, entitled PSIM (Photorefractive SIMulation), only accounts for the 
nonlinear photorefractive effect, not for SPS or other competing nonlinear effects. In 
so doing, I isolated the photorefractive effect from other effects and determined the 
necessary change in index of refraction across the beam path to account for beam 
fanning which is typically observed in the laboratory. The source code for PSIM is 
found in Appendix B. 

B. THEORY. 

As discussed above, the photorefractive effect may be solely responsible for the 
observance of beam fanning, which is in turn responsible for self-pumped phase 
conjugation [Ref. 2]. The change in index of refraction of a given crystal. An, may be 
given by the equation [Ref 14] 

An(z)= -(V 2 )n%frE(z). (3-1) 

where n is the index of refraction, r^^ is the effective Pockels coefficient and E(z) is the 
space charge electric field. 

Values for r^^p range generally in the picometers per volt range. Of seven crystals 
used by Klein and Schwartz [Ref 19], r^jp ranged from 4.2 x 10*^^ to 12.0 x 10'*^ 
meters per volt. Although these values are very small, the space charge electric field 
may be very high, on the order of 10^ volts per meter [Ref 5: p. 418]. Using these 
values it becomes obvious that the change in index of refraction due to the 
photorefractive effect in barium titanate may be as high as 10"^ . Though not an exact 
number for any given crystal, this may be accepted as close to an upper bound for An 
(denoted here as An^^^^^) for most of the crystals used. Gunter [Ref 11: p. 226] gives 
An^^ax ^ limited by the possible space charge field. 
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From the above simple calculations it may be assumed that for the 
photorefractive effect to be solely responsible for self-pumped phase conjugation, the 
dimensions of any given crv-stal used must be such that a change in index of refraction 
on the order of lO"'*^ is sufficient to cause the observed beam fanning. However, as 
previously mentioned, the gradient of the electric field and hence the index gradient is 
the important quantity relating to the photorefractive effect. Therefore, the beam size 
relative to the length of the crystal along the direction of propagation and the intensity 
profile must also play an important role in SPPC. These factors, in addition to the 
incident angle, were eventually incorporated into the simulation described here. The 
resulting simulation was used to determine the minimum change in index of refraction 
necessary to account for observed beam fanning in any given configuration. 

C. THE SIMULATION 
1. Program Logic 

The simulation is written in HP basic and designed to run on a Hewlett 
Padcard 300 model computer with at least two megabytes of random access memory. 
This two dimensional model takes into account only the change in index of refraction 
across the incident beam path due to the photorefractive effect. The output is a 
graphic representation of the eventual steady state beam geometry. 

The model uses a brute force approach to calculate the observed effects which, 
although time consuming and memory intensive, uses a simple algorithm. The 
program uses only Snell's law and Huygen's principle to calculate the internal beam 
pattern. Input into the program is the range of the change of index of refraction 
across the input beam. An, the crystal and beam size (in arbitrary units), the point of 
entrance of the beam into the crystal, and the angle of incidence of the incident beam 
relative to the normal. The program computes the mean index of refraction between 
two input extremes which then corresponds to the normal index of refraction (n ~ 2.4). 
The important factor in the program is the difference between the two extremes, 
"^^max • Logically it is this entity which drives the photorefractive effect. Note that 
che simulation is not concerned with how An^^,. is created. 

Once the necessary parameters are entered, the internal angle of the incident 
beam is calculated and each point on the wavefront is then assigned an index of 
refraction value that it can impart to the ciystal at its location. These points are then 
grouped into pieces of the wavefront (one to five separate pieces are normally used). 



32 



The assigned index value is dependent on and an assigned profile which may be 

varied (more will be said about this profile later). The indices of each point of the 
wave front are then propagated in a straight line through the crystal, at which time the 
value of the index of each point in the wavefront is imparted to the corresponding 
point within the crystal. After completing this process the w'avefront is propagated 
along the same path, allowing each point to move at a speed proportional to the index 
at the point in the crystal in which it finds itself After a short time the process of 
wavefront movement is stopped and the direction of each- piece of the wavefront is 
calculated using Huygen's principle. The original indices of the points within any given 
piece are then propagated in that new direction, imparting those values along the new 
path of the piece of the wavefront and the process begins again. The result of many 
iterations of this algorithm is a graphic output of the steady state internal beam 
pattern due to the photorefractive effect. 

2. Limitations 

There are several limitations to the model that must be discussed prior to the 
introduction of results. The first, and major, limitation is that of necessity each point 
on the wavefront cannot be treated individually. The wavefront must be treated in 
sections to determine an appropriate direction for a section and this makes the model 
only an approximation to the actual continuous case. The choice of the number of 
sections is primarily dependent on the chosen index profile. A profile similar to that 
shown in Figure 3.1 may be treated as a single piece since the direction of each 
individual piece will be the same as the direction of the entire wavefront. By contrast 
compare the profile of Figure 3.2. This profile appears to lend itself to the three pieces 
shown; however, in using only three pieces the effect around the extremum are 
minimized and it may be wise to consider the wave in five pieces. 

The difference in considering three as opposed to five pieces in the wavefront 
may be seen in Figures 3.3 and 3.4. In Figure 3.3 the index profile shown in Figure 3.2 
was used in three pieces. Note the distinctness of the three discrete sections. Figure 
3.4 shows the result of the same profile divided into five sections. It is still easy to 
note the five individual sections under consideration, however the output Is noticeably 
different. Ideally it would be desirable to have a very large number of pieces used in 
wavefront calculations. However, the larger the number of pieces, the larger the 
number of points in the wavefront, and consequently the larger the beam size. To keep 
the beam/crystal ratio in the proper perspective it is necessary to increase the size of 
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Figure 3.1 Linear index profile considered as one piece. 

the crystal whenever increasing the size of the beam, and the maximum cr>'stal size is 
limited by the available memoiy'. The maciiinc used for this research was limited to 
two megabytes of random access memory, uhich allowed for a maximum ciystal size of 
500 units. To keep appropriate beam/crystal ratios, the beam size used (generally) was 
between 20 and 40 units and the number of pieces in the wavefront ranged from one to 
eight, depending on the index profile. 

A second limiting factor in the use of this simulation is seen in cases where 
there is beam crossover. This happens when the index profile has a larger gradient on 
the side away from the c-axis than it does on the c-axis side. An example of this type 
of profile can be seen in Figure 3.3. In this case the beam intensity profile is 
considered to be Gaussian and the more incense section of the beam (i.e. the portion 
nearest the center) dominates over the weaker section. This ignores any consideration 
of two beam coupling which would actually occur. 

The final major limiting factor is the inability of the model to fully account for 
a piece of the wavefront which breaks away from the main portion of the incident 
beam. A case where this happens is seen in Figure 3.3 . In this case the index of 
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Beam Cantor 



Figure 3.2 Nonlinear profile with possible divisions. 

refraction gradient across the piece is the same as when it was part of the entire beam. 
This may be a proper assumption if the beam was continuous during fanning; however, 
as can be seen in Figure 2.4 and many others presented here, there are often portions 
of the beam which leave the major portion of the incident beam entirely. Under the 
condition of a complete disconnection it is possible that the total index gradient should 
be based on and not on the index gradient as it was calculated when ’the 

wavefront piece was within the incident beam; this was not programmed iiito the 
model. There are presently no data available indicating the effect of a detached piece 
of a beam in contrast to the continuous case. This problem is circumvented when a 
profile such as that shown in Figure 3.1 is used so that only one wavefront piece is 
considered. 

D. RESULTS 

As noted above, the output of the simulation is dependent on the index profile, 
the beam/crystal ratio, the angle of incidence and the position on the cr>'stal where the 
incident beam enters. The criteria for success (of possible SPPC) of a given 
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Figure 3.3 PSIM output for profile shown in Figure 3.2. 

configuration was taken to be the ability of a portion of the beam to reflect off the 
crj’stal corner in such a way as to rcintersect itself This would allov^' for, the 
production of auxiliarj’ beams which are a prerequisite to SPPC. Therefore, a result 
such as that seen in Figure 3.5 would be considered capable of SPPC, and one such as 
that seen in Figure 3.6 would not. 

Since the purpose of this investigation is to model beam fanning, the most 
important parameter is the inde.x. profile. To concentrate on this parameter, a 
beam, 'crystal ratio and angle of incidence that matched one of the configurations seen 
in the laboratory were chosen. The configuration chosen was an often used 
configuration with a beam/crystal ratio of 0.064 and an incident angle of 40° . The 
incident beam entered the crystal at a point way across the face of the crystal. This 
corresponds to a cry’stal observed in the laboratory of 5.0 mm on a side with a beam of 
size 0.32 mm, incident 3 mm from the front side of the crystal. 
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Figure 3.4 PSIM output for profile shown in Figure 3.2 
but broken into five separate pieces. 

Modeling was run with three index profiles. The rationale behind the use of the 
profiles will be explained in Section E of this chapter. In each case the values of 
AUmax which caused successful SPPC due to beam fanning are presented. 

Following the theory of Feinberg [Ref. 2], the first index profile used was that 
shown in Figure 2.5. Successful SPPC from beam fanning was seen (Figure 3.8), 
however, note the beam fanning in the direction opposite the c-axis; this is never 
observed in practice. 

The second beam profile used is shown in Figure 3.1. A typical successful run is 
shown in Figure 3.9. Although this approximates some observations of beam fanning, 
note the lack of a continued incident beam as one often (but not always) observes. 
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Figure 3.5 Beam fanning resulting in successful SPPC. 

The final profile used is shown in Figure 3.7. SPPC due to beam fanning was 
exhibited as is seen in Figure 3.10. It was this configuration which exhibited results 
closest to those observed in the laboratory. 

E. DISCUSSION 

The first index profile is taken directly from Feinberg [Ref. 2]. To arrive at a 
profile such as this it is necessary to assume that charge migration in the 
photorefractor is symmetric, resulting in a net charge distribution such as shown in 
Figure 2.5. The beam pattern predicted by the simulation is shown in Figure 3.S . The 
major problem with this profile is the fanning of the beam in the direction opposite the 
c-axis due to the tail of the index profile. In his original work Feinberg only considers 
the central region since it is the most intense part of the beam, ignoring the tails of the 
profile. As can be seen from the model, the tails of the profile cannot be totally 
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Figure 3.6 Beam fanning not resulting in successful SPPC. 



ignored. To arrive at the observed effect with the index profile of Figure 2.5 it is 
necessaiy to assume that barium titanate exhibits asymmetric self-focusing (i.e. self- 
focusing only in the direction of the c-axis). If asymmetric self-focusing was present, 
then the (less intense) light subject to the profile of the tails would be unable to fan 
against the c-axis. To date there has been no reported evidence of self-focusing of any 
kind in BaTiO^; however, this does not preclude the possibility since it would only 
affect a weak, portion of the beam and may not be obvious to the casual observer. 
Note that the fanning against the axis on the left-hand side of the beam would be 
abated by two beam coupling and might not be readily discerned. 

An alternative explanation is that the crystal response to the beam is not 
characterized by the curves shown in Figure 2.5. Instead, all available charge may be 
totally excluded from all portions of the beam, resulting in the profiles shown in Figure 
3.11. 
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Figure 3.7 Index profile for simulation of beam fanning in BaTiO^. 

In an effort to approximate an index profile which would exhibit the necessary 
behavior, the profile shown in Figure 3.1 was used. This beam profile was arrived at 
by making the important assumption that there is a total exclusion of charge carriers 
within the illuminated section of the crystal. This is an easy profile to work with since 
the outcome is independent of the number of pieces considered in the wavefront, as 
previously mentioned. The result of the simulation with this profile is shown in Figure 
3.9. However, there are two major problems with this profile. The first one is the lack 
of extension of the incident beam as mentioned above. Note that in mosi of the 
photographs presented in Appendix A there is an extension of the incident beam. This 
is almost universally, though not exclusively, true (see Figure A. 5 for example). 
Second, it is unrealistic in its sharp transitions. Both of these problems are solved wuh 
the introduction of the third index profile. 

The final index profile (Figure 3.7) resulted in a beam geometry which exhibited 
all of the requisite properties. This profile wa.s arrived at by the rounding of the 
discontinuities of the prior profile; specifically, a sine function was used to arrive at the 
configuration shown in Figure 3.7. The results of the simulation are presented in 
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Figure 3.8 PSIM output using the index profile of Figure 2.5. 

Figure 3.10. To arrive at this profile it would be necessary for diffusion of the charge 
carriers almost completely out of the incident beam path. The eventual diffusion of 
charge carriers completely out of the path may actually be expected since any light of 
intensity above ambient will cause any charge carriers to eventually be totally 
liberated.'^ The graphs of Figure 3.11 show the necessary charge, electric field and 
intensity profile of the incident beam necessary to arrive at an index profile similar to 
that used in this instance. It is assumed that the continuation beam on the left would 
be absorbed into the fanned portion by two beam coupling, eliminating the split 
appearance of the continued beam. 



■^Feinberg shows however, that the electric field may be proportional to the 
gradient of the incident beam intensity [Ref 27]. 



41 



To arrive at the configuration of Figure 3.9 a value of = 0.03 was used. 

In Figure 3.10 An^^^,^ = 0.02. Regardless of the profile used, these values of Arij^^x 
are characteristic of the value necessary to account for the observed beam fanning in 
BaTi03. For successful SPPC with the configuration chosen, An ^^^ must be greater 
than 10*^. Reducing the beam size will obviously reduce the necessary value of 
Anmax> values on the order of 10'^ to 10'^ are typical. 

Thus, the simulation has predicted the necessary change in index of refraction 
across the beam path, and the necessary profile of that change. In Chaper V the 
results of the necessary value for Anj^^j^ will be compared to those possible for the 
crystal. Also, the physical implications of the necessary profile of the index of 
refraction within the beam path, as detennined by the simulation, will be discussed. 
The comparison of the information derived from the simulation to that experimentally 
derived from the crystal, will show that the necessary conditions cannot be met in 
BaTiOj. 
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Figure 3.9 PSIM output using the index profile of Figure 3.1. 
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Figure 3.10 



PSIM output for index profile shown in Figure 3.7. 
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Figure 3.1 1 N'ecessary charge distribution and electnc field 
to produce the index proiile of Figure 3.7. 
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IV. EXPERIMENTS 



A. INTRODUCTION 

In order to investigate beam fanning in BaTiO^, an experiment was designed to 
measure the influence of a change in the incident beam diameter on the fanning 
process. If fanning increases with a decrease in beam diameter, it could be concluded 
that beam farming is a direct result of the change in index of refraction within the 
beam path, as predicted in Equation 2.25. If, however, the fanning decreases with a 
decrease in incident beam diameter, it could be concluded that beam fanning is actually 
forward-scattering due (probably) to SPS, since, as shown in Equation 2.28, the gain 
per unit length will decrease with a decrease in spot size. 

B. EXPERIMENTAL APPARATUS 

The apparatus used for this work is shown in Figure 4.1. The laser was a Spectra 
Physics Ar”^ variable power (20 mW maximum), multimode laser, model #162A-07. 
During the experimentation the laser was operated at all of the visible lines and was 
not optically isolated from the experiment. Upon leaving the laser, the beam was 
directed into a polarization rotator in order to rotate the beam from the ordinary to 
the extraordinary polarization in the crystal. To ensure that only extraordinarily 
polarized light was used, the beam was then directed through a Glan-Thompson 
polarizer. Following the polarizer was a beam sphtter, to enable the observation of 
any phase-conjugate reflection, and a variable aperture of the type common in 
photographic lenses. This aperture acted as a spatial filter, reducing the amount of 
scattered light reaching the crystal. The beam was then focused onto the crystal by 
uncoated lenses of either 34 cm or 16 cm focal length. 

The crystal was one of four single ciy'stals of BaTiO^ manufactured by Sanders 
Associates and poled into a single domain. Two of the crystals (on loan from the Los 
.-vlamos National Laboratories) measured approximately 5x5x5 mm and were 
transparent with a slight yellow tint. Two of the crystals (belonging to the Naval 
Postgraduate School) measured approximately 5 x 5 x 2.5 mm and had no noticable 
tint to them. The crystal was mounted on a rotatable platform capable of being 
translated 12 cm along the axis of the optical bench on which all of the components 
were mounted. Following the crystal, and connected to the mounting table, was a 
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knife edge. The knife edge was coated with black tape to ensure that a resonating 
cavity was not formed \\ith the crystal. Behind and to the side of the knife edge, also 
connected to the translation table, was a silicon photodiode connected to an 
electrometer measuring the short circuit current. 

C. EXPERIMENTAL PROCEDURE 

To measure the amount of beam fanning quantitatively,- the knife edge was 
aligned in such a way that the beam transmitted through the crystal was blocked from 
the detector. To ensure that the blockage of the transmitted beam was sufficient, the 
table was translated rapidly along its entire range of motion and the detector output 
was monitored to ensure that there was little or no change in the recorded intensity. 
The vibration inherent in the movement and the changing beam diameter of the 
incident beam ensured the lack of significant beam fanning, so that a lack of change in 
the recorded intensity indicated that the transmitted beam was completely blocked 
from the detector for all observed beam diameters. 

Once the crystal, knife edge and detector were in position, the beam was blocked 
and the crystal was flooded with light from a 75 watt incandescent lamp located 
approximately 12 cm away. After approximately one minute the light was turned off 
and the laser beam was unblocked. The intensity of the light reaching the detector was 
then monitored and the peak and steady-state intensities were recorded. Once the 
intensity had reached a steady-state value the beam was again blocked, the crystal was 
translated and flooded with white light, and the procedure was repeated. 

The flooding of the crystal with white light effectively homogenized the charge 
distribution within the crystal and ensured that one reading was not affected by the 
previous one. Special care was taken to ensure that the beam was parallel to the 
crystal table and that all components were mounted on the optical axis of the optical 
bench. The precision translation mechanism on which the crystal table was mounted 
ensured that the beam was incident on the crystal with the same center for all 
measurements. These precautions effectively compensated for any surface 
inconsistencies on die cn-'stal face except at very small laser beam diameters. 

Both peak and steady-state values of the fanned beam intensity were recorded; 
however, only the peak values were used for final analysis. This decision was made 
because, although the steady state values generally correlated with the peak values, the 
steady state values could be affected by numerous variables. One obvious problem 
with steady-state readings is the possibility of phase conjugation, thus reducing the 
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amount of fanning recorded in the steady state. Also, due to the orientation of the 
detector, it was possible (even likely) that the beam would fan across the detector, 
eventually coming to rest past the area it can record; sometimes the fanning was 
observed to exit the face perpendicular to the face the detector was monitoring. To 
reduce all of these effects, the incident beam was originally allowed to enter at only 
small angles to the c-axis ( ~ 2® - 5°) and, although alignment was tedious, consistent ^ 
results were obtained for various crystals, lens focal lengths, laser lines and incident 
angles. Eventually larger incident angles were used (~20° - 40°) with similarly 
consistent results. 

To ensure that the observed effect is not influenced by scattering initiated by 
small (e.g. molecular size) crystal inhomogeneities, one configuration was chosen in 
which there was moderate fanning and the fanned intensity was recorded for each of 
the available visible lines of the laser. The result was a mean ratio of incident intensity 
to fanned intensity that was constant across all available wavelengths, thus indicating 
no resonances and no dependence of scattering on wavelength at the available laser 
lines. 

After the data w'ere collected, the beam diameter at each point was measured 
with an Aeronca Electronics Laser Blade and the direct laser intensity was recorded. 
Here, the beam diameter is defined as twice the distance from the center of the beam to 
the point of 10% intensity. For analysis purposes the beam diameter could be divided 
by the cosine of the angle of incidence so that the beam diameter would be that 
actually along the y-axis of the crystal. 

Measurements were also made by placing the detector above the crystal as the 
incident beam diameter was changed. This gave an indication of the total amount of 
scattering, not just asymmetric scattering. In this case, the front surface of the crystal 
was shielded from the detector with black tape to keep specular reflection off of surface 
irregularities from influencing the measurements. Transmission measurements were 
also made by replacing the knife edge with an aperture, thus allowing only the 
transmitted portion of the incident beam to reach the detector. 

D. RESULTS 

Typical results of the fanning measurements are shown in Figure 4.2. In this 
case the cr>'stal was at an angle of < 5° so the beam diameter need not be adjusted for 
cos6. All intensity recordings have been normalized for comparison purposes. The 
shape of the curves is insensitive to the laser line used, beam angle, choice of crystal. 
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and the divergence or convergence of the incident beam (i.e. whether the crystal is 
placed in front of or behind the focal point of the focusing lens). 

Figure 4.3 is a plot of the intensity seen from above the crystal for various beam 
diameters for converging light. There are two important things to note here. First, 
there is less scattering for smaller beam diameters and second, the sharp discontinuity 
at about three micron beam diameter. This discontinuity is noticable as the crystal 
goes through the focal point of the focusing lens. Figure 4.4 is a plot of the scattering 
as viewed from above the crystal for diverging light. Note the initial decrease in the 
scattering within the first four microns. This is only apparent with diverging light and 
is seen in all crystals I used. Note also that there is a much broader dynamic range in 
the intensity for diverging light than for converging light. 

Although all of the scattering data are generally insensitive to which crystal is 
used (i.e. the general shape of all the curves for a given situation is the same for all 
crystals), the same cannot be said for the transmitted beam. Figures 4.5 and 4.6 show 
the transmitted intensity for all four crystals plotted together. Crystal #2 is the only 
one which does not exhibit SPPC and is the only crystal in which the transmitted beam 
intensity increases with beam diameter as a general trend. However, note that all of 
the curves exhibit some inflection point around six micron diameter for diverging light. 
This inflection is not a consistent maximum or minimum; it depends on the crystal. 
Presently I cannot explain the cause of these dissimilar curves. Generally it may be 
noted that the transmitted intensities map the same characteristics as the scattered 
intensity, but inverted, as would be expected. All of these various features will be 
addressed later, but the important aspect of all of the data is obvious; increased beam 
diameter results in increased beam fanning. 

Although these data alone are a strong indication of the lack of index change 
within the beam path due to the photorefractive effect, there exists even stronger 
evidence. By removing the knife edge and detector from behind the crystal, the image 
of the transmitted spot could be observed on a screen approximately three meters 
away. Figure 4.7 shows the image of the spot, which was projected directly onto 
photographic film at different times after the beam was allowed to enter the crystal, 
and prior to the onset of phase conjugation. Figure 4.7a is an image of the beam 
immediately after turn on; notice the spot is generally undistorted. Figure 4.7b is a 
photograph betw'een turn on and steady-state, and Figure 4.7c is a photograph at 
steady-state. By viewing these images it is clear that the energy of the beam is first 
extracted from the outside of the beam and not from the inside. 
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In addition to the above measurements, the index of refraction of the crystal was 
estimated by measuring Brewester's angle. At X = 514 nm the published value for the 
extraordinary index of refraction is n^ = 2.424 [Ref. 22], however, my measurements 
consistently resulted in n^ = 2.26 ±0.02. This large discrepancy is still unexplained. It 
is also interesting to note that, although there was a distinct minimum at Brewester's 
angle, the reflection did not disappear entirely. The minimum was so slight that it was 
quite difficult to find, and eventually a photodetector had to be used to accurately 
locate it. It is possible that both of these observations may be explained by the 
presence of impurities on the surface of the crystal. Specifically, an oil film from the 
hands of the experimenter may be at fault. Additionally, there is the possibility that 
heating at the surface of the crystal due to the high energy density of the incident beam 
may be responsible. Since measurements by others [Ref 30] of the index of refraction 
of BaTiOj using Brewester's angle have resulted in indices consistent with the reported 
values, either one or both of these explanations may account for the observed 
anomalies. 
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Figure 4.1 Experimental Apparatus. 
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Figure 4.2 Dependence of fanning intensity on incident beam diameter. 
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Figure 4.3 Typical dependence of scattered intensity on beam diameter 
as measured from above for converging light. 
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Figure 4.4 Typical dependence of scattered intensity on beam diameter 
as measured from above for diverging light. 
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Figure 4.5 Transmitted intensity as a function of beam diameter for diverging light. 
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Figure 4.6 Transmitted intensity as a function of beam diameter for converging light. 
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Figure 4.7a 



Photograph of projected spot immediately upon turning on the incident beam. 
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Figure 4.7b Photograph of projected spot after the onset of beam fannin 
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Figure 4.7c Photograph of projected spot in the steady state. 
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V. ANALYSIS 



A. DISCUSSION OF RESULTS 

The simulation in Chapter III and the experiments described in Chapter IV 
clearly show that asymmetric self-defocusing cannot be due solely to a change in index 
of refraction along the beam path within the crystal. Even with the changes of the 
profile for the change in index of refraction proposed in Chapter III, the necessary 
value of is too large according to the published maximum values for BaTiO^. 

Inspection of the necessary values of An^^j^ from the simulation shows that values on 
the order of 10'^ are necessary to account for the observed fanning within the crystal; 
however, previously calculated values place an upper bound of ~ lO"^, and a value of 
~ 10'^ is probably more reasonable^ [Ref 11]. 

This discrepancy may be accounted for only by assuming that the index gradient 
is limited to a very small portion (in some cases as small as the size of the input 
beam) of the center of the beam, or by assuming that the observed effect is a result of 
another process. It is impossible to justify the first assumption and stiU achieve the 
necessary index profile, since charge carriers must be liberated from the entire region 
where the brightness of the incident beam exceeds ambient. Therefore, by Laplace's 
equation, the gradient of the electric field within the beam path must be uniform. This 
leaves only the possibility that another process is responsible for beam fanning. 

The results of the experiments of Chapter IV reinforce the conclusions made 
from the simulation portion. Figure 4.2 shows that beam fanning decreases with a 
decrease in beam diameter. From Equation 2.25 it is clear that this behavior cannot be 
explained by a change in the index of refraction within the beam path due to the 
photorefractive effect. If this type of change were responsible for beam fanning, the 
the amount of fanning would have to increase with a decrease in beam diameter. 
Additionally, as shown in Figure 4.7, beam fanning begins on the periphery and the 
light IS scattered progressing from the outside of the beam toward the center. Noting 
the index profile as shown in Figure 2.2, this is impossible if the fanning is due to an 
asymmetric change in index of refraction across the beam path, as was originally 
predicted. 

^Feinberg does, however, calculate a possible maximum value on the order of 
10‘^ [Ref 27]. 
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One may argue that all of the observed qualities of beam farming may be 
explained by assuming "small-scale" self-focusing of the incident beam within the 
crystal. That is, the incident beam may experience a breakup into spots which are of 
diameters small enough to create index gradients of the magnitude shown necessary by 
the simulation in Chapter III. This type of self-focusing has been known for over two 
decades, and is common in materials whose dielectric constant increases with field 
intensity [Ref 31]. If small-scale self-focusing occurred in BaTiO^, it could explain 
many of the effects that are observed during SPPC (e.g. the presence of small filaments, 
beam farming, etc.). However, there are serious problems with a small-scale self- 
focusing theory. First, the breakup of optical beams due to self-focusing has been 
shown to be spatially periodic, with the periodicity being dependent upon intensity 
[Ref. 32]. As is demonstrated in Figure 4.7, there is no indication of spatial periodicity 
in the beam as it is transmitted through a BaTiO^ crystal. Also, even in good self- 
focusing mediums at very high beam intensities ('*^ 10^ W/cm^), it is necessary for the 
beam to have a long path length (on the order of tens of centimeters) within the 
medium before self-focusing is observed [Ref. 33]. Even with a very large change in the 
dielectric constant, it is unlikely that self-focusing could occur within a medium of only 
5 mm in length. Finally, there is no evidence in the literature that the dielectric 
constant of BaTiO^ changes significantly with the intensity of the incident light 
(although it does change with an intensity gradient, as discussed in Chapter If). It is 
interesting to note however, that the dielectric constant of BaTiO^ does change with 
temperature [Ref 25], and this may cause self-focusing at very high powers and very 
small beam diameters, given a long enough path within the medium. 

From the evidence presented here it follows that asymmetric self-defocusing in 
BaTi 03 cannot be due to an asymmetric change in index of refraction along the path 
of the incident beam due to the photorefractive effect. Attributing this phenomenon to 
some other process however, is not so easy. 

B. SCATTERING 

As noted in Chapter II. the only other theory that is currently available to 
explain asymmetric self-defocusing in BaTiO^, besides the theory just discussed, is the 
theory involving scattering (SPS) [Ref 24]. Although I have not conclusively 
demonstrated that beam fanning in BaTiOj is due to SPS, there are a number of 
observations which indicate that this type of scattering maybe the origin of the effect. 
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As may be seen in Equation 2.28, if near-forward scattering is the origin of beam 
fanning, then there should be a decrease in fanning with a decrease in incident beam 
diameter, an observation borne out in our experiments. Also, the diversion of energy 
from the incident beam could progress from the periphery toward the center, since 
noise photons originating from outside the beam must interact with the first portion of 
the beam they encounter (Figure 4.7). 

Another interesting and important discovery is that although the shape of the 
fanned intensity curves is similar regardless of convergence or divergence of the 
incident beam (in front of or behind the focal point of the focusing lens), the absolute 
intensity of the scattered light is always observed to be greater for the diverging case. 
This is not perplexing when considering near-forward scattering. One would expect 
more scattering from diverging light than converging because the size of the beam 
increases as it progresses through the crystal for diverging light, thus increasing the 
amount of scattering. Also it is noted that individual intensity readings are not 
consistent, even if there is no no movement of the apparatus between readings. This is 
expected since the scattering process begins with noise photons, and no two situations 
can be exactly alike, even for similar configurations. This inability to exactly reproduce 
any given measurement suggests that the scattering may begin with random index 
changes within the crystal, caused by the stochastic distribution of charge carriers. 

Finally, if stimulated scattering is responsible for beam fanning, then one would 
not expect continuous fanning, but fanning only along paths from which an incident 
photon interacted with the incident beam, thus creating a set of filaments emanating 
from the incident beam. A close examination of any of the photographs of Appendix 
A shows that this is indeed the case. Fanning is a set of discrete small beams all 
emanating from the incident beam. We must now discuss curved appearance of the 
dominant-scattered beam during self-pumped phase conjugation, a phenomenon not 
explained by SPS. 

A very close examination of many photographs of SPPC (e.g. Figure A. 5) reveals 
that beams which appear to curve often have a series of discontinuities within them. 
Generally, a filament, or set of filaments, leaves the incident beam along the direction 
of the c-axis. After a certain distance (usually about halfway to the crystal edge) these 
filaments enter a region where they become fuzzy, and often two become 
indistinguishable. The filaments exit this fuzzy region at an angle different than their 
entrance angle. This type of behavior is explainable by assuming that within this fuzzy 
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region there is stimulated scattering due to SPS. Forward scattered beams enter into a 
region where scattering occurs due to photons arriving from a different angle. The 
beams which continue on at a new angle are the forward scattered beams from this 
second scattering. It is therefore a series of this type of scattering discontinuities which 
gives the characteristic curved appearance of a filament during SPPC. However, this is 
not always the case. 

There exist configurations where this series of discontinuities does not appear to 
account for the curved filaments (i.e. Figure A. 6). In a majority of these cases a close 
examination reveals that a large filament originates at the incident beam or a series of 
filaments leave the incident beam very close together. The curved appearance comes 
from the coupling of energy to the filament, or portion of the filament, in the direction 
of the c-axis as the filaments progress through the crystal. Therefore the beam doesn't 
actually curve; it simply appears to do so because it is losing energy from one side and 
transferring it to the other. The process is actually a form of self-focusing due to two 
beam coupling. 

I have observed a very few instances where it appears that neither of these 
processes can fully account for the curved appearance of the fanned filaments. In 
these cases one would expect the effect to be due to true beam fanning from the 
change in index of refraction within the filament path. Although the simulation PSIM 
shows that it is impossible for this phenomenon to fully account for the observed 
fanning effect, it would be surprising not to find some evidence of a change in index of 
refraction on the order of 10'^ across a beam of diameter less than a few microns. 

C. SELF PUMPED PHASE CONJUGATION 

From the discussion of the previous section we may conclude that a change in 
index of refraction within the incident beam path is not responsible for asymmetric self 
defocusing in BaTiO^. It is probable that SPS is responsible for the effect. Therefore 
one may say that the true origin of SPPC in BaTiOj is scattering, probably SPS. The 
cause of the final phase conjugate reflection is not as clear. 

Phase conjugation by stimulated scattering in BaTiO^ has been demonstrated 
[Ref 23], and recently it has been shown how backscattering may dominate forward 
scattering in SPS [Ref 24]. Additionally, the investigations reported here indicate that 
scattering is responsible for the phenomenon of beam fanning. Therefore one might 
incorrectly conclude that phase conjugation in BaTiOj is due to stimulated scattering. 
With the onset of backward scattering in the crystal there begins a competition of the 
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scattered waves. As the modes begin to build, a mode showing high gain can usurp 
power from one showing less gain. In this way only the modes showing the highest 
gain eventually survive. It is this competition between configurations having high gain 
which is responsible for the amplitude oscillations often seen in a phase conjugate 
reflections (see [Ref. 34,35] or Appendix A). 

If a filament, arising from near forward scattering, has a configuration such that 
it can reflect off the comer of the crystal (i.e. a CAT comer) to create a situation 
where degenerate four wave mixing (DFWM) is possible in addition to SPS, then one 
would expect that filament to have a gain advantage. This then would be the filament 
to eventually win the competition. So it actually would be expected that both DFWM 
and SPS should contribute to SPPC. 

A review of the conditions of the only published report of phase conjugation due 
solely to SPS [Ref. 23] shows that special precautions were taken to ensure maximum 
coupling. Specifically, an index matching liquid was used which preserved as much of 
the phase conjugate return as possible. Additionally the crystal was surrounded by a 
glass cuvette (containing the index matching liquid) w'hich increased the amount of 
scattered light available for SPS to begin. Even with these advantages, a phase 
conjugate reflection of only 10% was obtained. This is quite small compared to the 
30% to 50% usually seen from SPPC in BaTiO^ in air. This indicates that DFWM 
from the internal comer reflection is the most significant contributor to the phase 
conjugate return observed in SPPC, as was originally proposed [Ref. 2]. 
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VI. CONCLUSIONS 



Experimental evidence does not support the theory that asymmetric self 
defocusing in BaTiO^ is a result of an asymmetric change in the index of refraction 
within the path of the incident beam. Instead, forward stimulated scattering, probably 
Stimulated Photo refractive Scattering (SPS), is found to be a more likely explanation. 

Simulation has shown that for an index change to be responsible for the 
phenomenon of beam fanning, the necessary change in the index of refraction across 
the beam must exceed the theoretical maximum, based on the maximum possible space 
charge field. A change in the index of refraction on the order of 10'^ is necessary to 
account for the observed effects, whereas the accepted maximum is on the order of 10'“^ 
to 10'^ 

Experimental results do not support a theory which requires a change in index of 
refraction across the incident beam as the cause of beam fanning. The magnitude of 
fanning is observed to increase with an increase in beam diameter. Also, observations 
of the manner in which the beam fans show that energy is scattered out of the beam 
path beginning with the periphery and not beginning in the center. Neither of these 
observations is consistent with the results expected if a purely photorefractive effect 
were responsible for asymmetric self defocusing. 

Close observations of the process of beam fanning in BaTiO^ indicate that beam 
fanning can be explained by a stimulated scattering theory. The presence of filaments 
rather than continuous fanning, scattering beginning with the periphery of the incident 
beam, and discontinuous filaments during phase conjugation are all e.xplicable when 
beam fanning is attributed to stimulated scattering. 

Future experiments in this area should include probing the "fuzzy"' regions of the 
observed filaments from above with a low power laser to search for the presence of 
additional index gratings. Also, the actual change in the index of refraction across the 
incident beam path witliin the ciystal remains to be determined experimentally. 
Additionally, further research into the size, type and spacing of scattering centers 
within the crystal is necessary to fully explain the observed difference in scattering of 
diverging and converging light. 
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APPENDIX A 

INTERNAL BEAM PRODUCTION 



1. INTRODUCTION 

The production of internal beams in BaTiO^ has, to my knowledge, never been 
systematically investigated. A search of the available literature indicates that the cause 
for this oversight may be that, in general, millimeter sized, collimated beams have been 
used. In aU available photographs of BaTiO^ used as a SPPCM, the beam has not 
been focused onto the crystal. If the beam is not focused, then the diameter of the 
beam upon entering the crystal is at best the same as when it leaves the laser. It 
appears that when a beam is not focused onto the crystal, auxiliary beams are 
eliminated or obscured by the strong input beam since its diameter is large compared 
to the dimensions of the crystal (sometimes being more than half the size of the 
crystal). However, if the input beam is focused, the beam diameter is small compared 
to the crv'stal and the immediate effects of the input beam stay localized within the 
crystal, thus allowing the formation and observation of auxiliary beams. There is also 
an immense increase in the energy density of the beam within the crystal upon 
focusing. 

It was using a focused beam that we furst discovered the ability of barium 
titanate to support internal auxiliary beam patterns. Extensive investigation by A. V. 
Nowak and myself at Los Alamos National Laboratories indicated that there is a 
plethora of internal beams which may be generated by the nonlinear effects of barium 
titanate during phase conjugation. It appears that phase conjugation is necessary to 
sustain these beams since no auxiliary beams were observed in the absence of phase 
conjugation. In many cases it seems that the beam is initiated by the reflection of the 
phase-conjugate beam off of the face of the crystal as it exits, but this is not always 
true. 

What follows in this chapter is an attempt to report, categorize, and explain 
many of the observed internal beam patterns and anomalous effects. It is important to 
note that often it is difficult or impossible to recreate some of the observed effects. It 
appears that not only position and beam profile are important, but also the short term 
history of the crystal. This means that possibly one effect may not be observable 
unless another has (or has not) immediately preceded it. 
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2. EXPERIMENTAL APPARATUS 

All investigations presented in this chapter were performed at Los Alamos 
National Laboratory on an air floated optical bench in free air that was maintained at 
a temperature of 22± 1°C. The apparatus used is shown schematically in Figure A.l. 

During all experimentation the light source was a Spectra Physics 2020-03 argon 
ion laser, prism-tuned to 514.5 nm and operated in a single longitudinal mode by 
means of an intercavity etalon. The beam at the laser had a diameter of 0.88 mm at 
the Vq of peak intensity. 

After leaving the laser, the beam was directed through a half-wave plate that 
served as a beam attenuator when necessary. Following the half-wave plate was an 
apparatus for isolation of the laser from the phase-conjugate beam. The apparatus 
used was a faraday rotator, manufactured by Optics For Research Inc., consisting of a 
25mm long rod of Hoya FR5 glass surrounded by a set of permanent magnets. The 
isolator was set for optimal rotation of the plane of polarization of the beam at 45® 
and placed between two polarizers with axes oriented at 45® to each other. This 
configuration gave a backward attenuation of approximately 500. Without the isolator 
the phase-conjugate beam entered the laser cavity and caused large fluctuations in 
pow’er output as well as mode hopping. The polarization of the beam upon leaving the 
isolator was such that it would enter the crystal extraordinarily polarized, which is 
required for the formation of a SPPCM. 

After passing through the isolator the beam was expanded by a factor of 6.25 
with an expanding telescope. This allowed the insertion of transparencies into the 
beam when desired, enabling the verification that an observed beam is a phase- 
conjugate beam. Also included in the telescope apparatus was a 0.05 mm pinhole that 
could be removed at will. This pinhole provided the ability to make observations with 
a Gaussian beam (pinhole inserted) or with what was eventually termed a mottled 
Gaussian beam (without the pinhole). The mottled Gaussian was so named because, 
although the beam without the pinhole was generally Gaussian, it had a mottled 
appearance. 

Followng the telescope, the beam was directed through a beamsplitter with 67% 
transmission at 514 nm. This beamsplitter not only provided the means for viewing the 
phase-conjugate signal as it returned along the incident path, but also provided a 
reference beam that was used to look for frequency shifts of the phase-conjugate beam. 
To look for a frequency shift the phase-conjugate return was reflected from the 
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Figure A. 1 Experimental Apparatus. 
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beamsplitter and redirected to interfere with a portion of the incident beam reflected 
from the front side of the beamsplitter. The resulting interference pattern was then 
observed for evidence of movement. 

Upon passage through the beamsplitter, the beam, with power of 75 ± 5 mW, 
was either modified by a lens, slit, knife-edge or a combination thereof, or simply 
collimated and directed toward the crystal. As the phase-conjugate beam returned 
along the incident path, the beamsplitter reflected a portion of it and directed it to a 
large area silicon PIN photodiode reverse biased at 22 volts and terminated with a 100 
ohm resistor. The output of the detector was sent to a digital millivolt meter, an 
oscilloscope and a chart recorder. When an image was placed in the beam the detector 
was removed and replaced with a white screen for viewing. When frequency shift 
information was desired, the detector was replaced with a mirror and the two beams 
were directed toward a white screen for observation of the resulting interference 
pattern. 

The crystal was any one of four nearly cubic single crystal of BaTiO^ measuring 
5 ±0.2 mm on a side. The sides of the crystals were polished, but it is important to 
note that the sides were not exactly the same size and therefore the sides were not 
exactly orthogonal. The cr>’stals were mounted on a platform that could be raised, 
lowered, rotated and tilted. Observations of internal beams were made visually from 
above and were recorded by film or video tape. Both the still and video cameras used a 
f = 35mm Cannon FL lens mounted on a single track bellows extended to provide a 
full frame view of the crystal. 

Many of the photographs presented here were taken from the video tape directly 
off of a television screen using a 6x7 cm format still camera. The desired frame of the 
video was digitized and placed on a hard disk for recall and photographing, the 
numbers displayed in the photographs are a result of the equipment used to do this 
final photography and are not pertinent to the investigation. When this method was 
used, it was possible to make negative images of the desired photographs, which 
reproduce better in print, therefore some of the photographs presented in this work are 
m the negative. There is no other significance to a negative versus positive print. 

In this and in the following sections any reference to phase-conjugate reflectivity 
is the percent return of the beam actually incident on the crystal. The amount of 
specular reflection off of the face of the crystal (which can vary [Ref 36]) is not 
subtracted out of the calculation. All reflectances have an experimental tolerance of 
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±3%. In the photographs presented, except where noted otherwise, the incident beam 
enters vertically from the top and is visible only upon entering the crystal. 

3. RESULTS AND ANALYSIS OF DATA 
a. Experimental variables 

The number of variables involved in SPPC is immense, however I have tried to 
concentrate on those w'hich have emerged as the most important. The most important 
factors in the production of various beams differ, but they- generally may be considered 
as dependent on beam profile or the geometry of the beam-crystal interface. 

Lenses, knife-edges, slits and pinholes were used to change the profile of the 
incident beam prior to the beam entering the crystal. The effect of these modifications 
of the incident beam is especially obvious when the beam profile has a large gradient, 
such as when a knife-edge is inserted into the beam, since this creates a very clear 
delineation in the lighted and unlighted areas of the crystal and introduces new wave 
vectors into the beam. The effect of the lens to crystal distance on phase-conjugate 
reflectivity has been studied by A. V. Nowak and is still not well understood [Ref 37]. 

The geometry of the beam-ciy’stal interface is probably the most important of 
all variables in the system. Changes in angle of less than one degree and lateral shifts 
in the beam of less than 0.5 mm can cause drastic differences in internal beams 
configurations. The following examples are presented to indicate the drastic effects 
observed with these type of changes. 

The configuration shown in Figure A. 2 is a typical experimental configuration 
that produced a phase-conjugate signal from the BaTiO^ crystal. The phase conjugate 
reflectivity was 33%. Note the fanned beams at the lower left corner of the crystal. 

Figure A. 3 demonstrates the extreme effects made possible when changing the 
beam profile by placing a knife-edge in the incident beam. Unless otherwise indicated, 
the knife-edge always intercepted the beam on the right side and was placed behind any 
lenses. Note the extreme curvature of the internal beam observed when blocking half 
the beam with a knife-edge. This also reduced the phase-conjugate reflection to zero. 
At other ingles and intersection points, ’■.he introduction of a knife-edge produced a 
phase-conjugate reflection where there was none before. 

A large effect is also seen when relocating the beam entrance position by a 
lateral shift in the beam-crystal intersection. As the beam is moved across the face of 
the crystal, the phase-conjugate reflectivity will change, beginning with no phase 
conjugation at one end and ending with the same result at the other. In between the 
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Figure A. 2 Normal configuration for BaTiOj exhibiting SPPC. 

two extremes, the phase-conjugate reflectivity may range as high as 60% or be 
nonexistent depending on other factors. Similarly, a change in the angle of the beam 
entering the crystal will have a drastic effect on the phase-conjugate reflectivity and 
beam pattern. It is generally unpredictable and small changes in angle (~4®) have 
been observed to have no effect, to destroy a pattern and eliminate all phase-conjugate 
reflectivity, or to increase the phase-conjugate reflectivity as much as 500%. 
b. Complex Auxiliary Beam Patterns 
— L The Total Internal Reflection ( TIR) Ring 

One of the most interesting effects discovered was the TIR ring. Figure 
A.4 shows the first example, to my knowledge, of a totally internally pumped ring 
oscillation within a BaTi 03 crystal. This effect was achieved using a knife-edge to 
obscure half of a mottled Gaussian beam that was then focused through the f= 30 cm 
lens, 29.5 cm from the face of the crystal. The beam was incident at 40° to the axis of 
the crystal and the knife-edge was placed on the side of the beam nearest the crystal. 
The effect was occasionally obtained while obscuring the other half of the beam, and 
was once obtained with a vertically aligned slit. A. V. Nowak once obtained the ring 



71 




Figure A. 3 Effect of placing a knife edge in the incident beam. 

with just the focusing lens, indicating that the TIR ring may be dependent more on the 
availability of light scattered at a certain angle than on the incident beam profile. 

In the case shown in Figure A.4 the phase-conjugate signal took about 30 
seconds to appear. The beam within the crystal slowly moved to the corner until the 
configuration of the beam arrived as shown, but without the ring. Then, slowly, the 
ring appeared over a three second interval. As the ring appeared, the phase-conjugate 
reflectivity dropped from 76% to 67%. A drop in phase-conjugate reflectivity was 
observed in all examples of ring formation for those cases where the reflectivity was 
measured. In some cases, several nearly overlapping rings were observed, as in Figure 
A.4. The ring in Figure A.4 was steady, fading occasionally (on the order of once 
every five minutes) but always returning. It was not overly sensitive to vibration, 
although large vibrations would destroy the ring. Rotation of the crystal destroyed the 
ring but if the rotation was under then the ring would reappear. Translation along 
the cry'stal on the order of 0.5 mm in either direction did not disturb the ring; further 
translation would destroy it. For other configurations, the ring was only preserved for 
one direction of crystal motion (see Section 4). 
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Figure A.4 The TIR Ring. 



The direction of circulation within the rings was determined by extracting 
some of the light out of the ring into a second cry'stal of BaTi 03 . After producing the 
ring in one crystal, a drop of index matching liquid (methylene iodide (CH 2 I 2 ) diluted 
1:1 with 2-propanol, n 1.6) was placed onto the surface of a second crystal and the 
second crystal was moved into contact with the first. Since the loss into the second 
crystal was too great to allow the continuation of the ring, the operation was recorded 
by the video camera and the frame of the instant was of contact was studied at a later 
time. Study of the video revealed that the light within the ring could propagate 
counterclockwise, or both ways simultaneously depending on the particular ring. 
However, in the majority of cases, only the counterclockwise direction was observed. 
Figure A. 5 clearly shows this counterclockwise propagation and is exemplary of the 
typical situation found. 

The ring within the crystal is probably produced by two-beam coupling 
between the main beam and scattered light that fmds a closed loop via internal 
reflections around the crystal. This somewhat lossy TIR ring takes power from the 
incident beam, which explains the drop in phase-conjugate reflectivity seen immediately 
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Figure A. 5 Determination of the direction of the TIR ring 
by extraction of light with another crystal. 

upon its appearance. Of interest also is the difficulty in finding the TIR ring; in some 
crx’stals the ring pattern could be found in minutes. In others it would take many 
hours. Nowak discovered that in many cases, however, formation of the ring was 
promoted by a slight forward or backward tilt of the crystal (out of the plane of the 
table). 

The effect of crystal motion on the TIR ring was also investigated by 
Nowak [Ref 38] and it is interesting to note that the TIR ring appears to seek a square 
shape. It was observed that the ring will go from rectangular to square upon 
translation but not from square to rectangular. Also discovered were configurations in 
which the ring exhibited total shape invariance under translation. 

3. The Diamond 

The first of the diamonds, the sharp diamond, is illustrated in Figure A. 6 It 
is considerably easier to obtain than the full ring and has been produced wth a variety 
of input beam profiles and orientations. This is thought to be a simple case of two 
beam coupling between the main beam and an internal reflection. Usually when the 
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Figure A. 6 The Diamond configuration. 

sharp diamond appears in the crystal the phase-conjugate signal increased considerably 
(up to 50 times what it is just prior to its appearance). The sharp diamond probably 
begins as a specular reflection of the phase-conjugate signal, which reflects back to the 
incident beam by internal reflection. Since a second point for DFWM is created, one 
would expect the observed increase in the phase-conjugate signal upon the appearance 
of the diamond. 

A fascinating and unexplained variation of the sharp diamond is the diffuse 
diamond. The diffuse diamond, shown in Figure A. 7, begins in a configuration similar 
to the diamond shown in Figure A. 6, but as the crystal is laterally moved, a portion of 
the incident beam composing the sharp diamond stays to create a diamond in the right 
half of the crystal. The diffuse diamond is very weak compared to the other effects 
discussed, but it is easily visible, especially during its creation. This etTect is similar to 
the motion-invariant effects which will be discussed later. 

5. Internal Specular Reflection 

Often the specular reflection of the incident beam within the crystal 
behaves in unexpected ways. Figure A. 8 shows a case where there is visible specular 
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Figure A. 7 The Diffuse Diamond configuration. 

reflection of the input beam within the crystal creating a diamond but there is no 
reflection of the secondary beam off of the rear face of the crystal. Absorption of the 
reflected beam by two-beam coupling with the incident beam within the first few 
micrometers of the crystal edge is probably responsible for the lack of specular 
reflection in this and other similar cases. There is a similar theory which accounts for 
the lack of specular reflection of the phase-conjugate beam off of surfaces such as the 
input face of the crystal. Pepper (Ref 36] proposes that the lack of a reflection of the 
phase-conjugate beam is due to destructive interference between the incident beam 
reflection, and a beam created by the phase conjugation of the reflection of the original 
phase-conjugate beam. This has been borne out in experiment. However, as Figure 
A. 9 demonstrates, at certain angles the specular reflection is sufficiently intense to 
make many additional internal reflections. The boundarv' conditions necessary for the 
lack of a reflection have yet to be fully determined. 

4. Diagonals 

Figure A. 10 illustrates the diagonals. The diagonals faded in and out as the 
phase-conjugate reflectivity oscillated from 12% to 21% with a period on the order of 
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Figure A. 8 Example of lack of reflection of the incident beam 

within the cr>’stal. 

5 seconds. The reflectivity was strongest w'hen the diagonals were brightest. When the 
angle of incidence was increased to 56® a similar configuration appeared, but w'ith the 
phase-conjugate reflectivity oscillating from one to 15% and specular reflection off of 
the back face fading as the diagonals became brighter. At another angle of incidence 
in this configuration ( 62®) it was noted that increasing the power of the input beam 
could eliminate all oscillations of the phase-conjugate reflection while eliminating the 
diagonals. This configuration has also been seen without any measurable phase- 
conjugate reflection. 

c. Temporal Effects and Frequency Shifts 

Up to this point this work has considered only examples of the many internal 
beam patterns that can be obtained by van.’ing the beam profile and orientation. 
Observations of oscillations in the pnase-conjugate amplitude and of frequency shifts 
between the phase conjugate signal and the incident beam will now be addressed. 
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Figure A. 9 Example of multiple reflections within the crystal. 

/. Amplitude Oscillations 

Frequently, amplitude oscillations were observed in the phase-conjugate 
signal of many of the previously discussed examples. These include Figures A. 2, A. 4, 
A. 6, A. 8. A. 9. The oscillations, probably similar to those reported by Valley and 
Dunning in an external resonating cavity [Ref 39], were most obvious when a mottled 
beam was focused into the sample, but they also appeared in some configurations with 
all beam profiles. The observed oscillations often appeared to be chaotic, and much 
work has been reported recently concerning the observation of chaotic oscillations in 
BaTi 03 [Ref 34,35]. However, some of the observ'ed oscillations were very nearly 
periodic (periodic within 10% of some mean period) with periods ranging from 0.01 to 
5 seconds. 

Figure A. II shows the amplitude of the phase-conjugate signal obtained 
with a focused Gaussian beam incident on the center of the crj'stal 40° to the c-axis. 
These oscillations appeared immediately upon the initiation of phase conjugation. 
These w'ere maintained as showm for a period of about one minute and then abruptly 
changed to the oscillations shown Figure A. 12. Within a few minutes the original 
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Figure A. 10 The Diagonals configuration. 



oscillations reappeared. Fast oscillations (~10 Hz) amplitude modulated by slower 
oscillations ('^0.5 Hz), and oscillations with faster beats superimposed on them were 
also observed. Nowak has shown that the frequency of these amplitude oscillations is 
proportional to the intensity of the incident beam [Ref 37] and this has been confirmed 
by others (Ref 30). 

These amplitude oscillations of the phase-conjugate signal are probably a 
result of frequency shifts of those auxiliary beams that account for the formation of the 
grating. Frequency shifts of the auxiliary beams were proposed by Feinberg and 
Bacher (Ref 40] to explain observed frequency sweeping in the output of a ring 
resonator using BaTi 03 and in an unisolated Ar"^ laser producing phase conjugation 
in BaTiOj (Ref 41]. Thus it is expected that the phase-conjugate reflectivity would 
diminish as the frequency diiTerence between the incident and reflected beams becomes 
larger, since the four wave mixing signal is reduced for large frequency differences as 
shown by MacDonald and Feinberg (Ref 42]. The nature of and reasons for this 
frequency shifting within the crystal has yet to be explained. 
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Figure A. II Example of quasi-periodic amplitude oscillations 

in BaTiOj. 



2. Frequency Sweeping 

To determine the difference in frequency of the phase conjugate signal, a 
Michelson interferometer arrangement was used. A portion of the incident beam 
interfered with the phase-conjugate signal, as described above, and the resultant 
interference pattern was observed for signs of motion. Normally there was no 
movement of the interference fringes; however, in the case of some configurations 
involving amplitude oscillations (e.g. Figures A. 6 and A. 10) there was some evidence 
of a frequency difference between the two beams. 

Usually as the phase-conjugate signal began to rise, the interference fringes 
were steady. However, as the signal began to decline the fringes began to move, and 
the more the signal declined the faster the fringes moved, until eventually they were 
indistinguishable. This sweeping process repeated itself with each oscillation of the 
phase-conjugate reflectivity. The direction of the movement of the interference fringes 
(the direction of frequency shift from the frequency of the incident beam) changed 
arbitrarily, implying that positive and negative frequency shifts were possible. This 
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Figure A. 12 Amplitude oscillations which progressed from 
those of Figure A.l 1. 

correlation of amplitude and frequency shift of the phase-conjugate beam may have 
been explained by MacDonald and Feinberg when they showed that the four-wave 
mixing signal falls off smoothly as the frequency difference between the reflected and 
incident beam becomes larger than ~ 1 Hz. Although the mechanism for this 
frequency shifting is still under investigation, Lam [Ref. 3] has proposed that the bulk 
photovoltaic effect gives rise to an intensity-dependent frequency shift. This implies 
that the frequency is dependent on the amplitude instead of vice versa. Another 
explanation may be that there is consistent movement of the grating when the reflected 
pumpwave forms a loop (see Figure A. 2 for example). In this case each circuit around 
the loop will Doppler-shift the pumpwave and the observed effect may be expected. 

There was only one observed instance ivhere the frequency difference 
betw'een the two beams was non-zero and constant (no photograph available). The 
phase-conjugate return was frequency-shifted from the input beam, but there was no 
sweeping. That is, there was obvious movement of the interference fringes at a 
constant velocity. This effect resulted from a mottled Gaussian beam incident on the 
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crystal at an angle of 35° to the c-axis. However this is an anomalous effect that has 
yet to be repeated or explained. 

i. Influence of Returning the transmitted Beam 

This effect is similar to that reported by Feinberg [Ref. 27] and was 
discovered while trying to stabilize the amplitude oscillations of the phase-conjugate 
signal discussed above. The effect was seen using an arrangement that produced 
nonperiodic amplitude oscillations of the phase-conjugate reflection ranging from zero 
to 46%. In this configuration a 100% reflectance mirror was placed behind the crystal 
directing the transmitted beam back through the crystal, but not directly along the 
incident path. All oscillations stopped, leaving the phase-conjugate reflection steady at 
10 %. 

The explanation for this effect is probably found in the competition 
between beam patterns. In any instance, the observed beam pattern is the one which 
exhibits the maximum gain. When oscillation occurs between two competing beam 
patterns it is apparent that, for some reason, the gain of one is falling above and then 
below the possible gain of another. The input of a beam of light from any source with 
brightness on the order of the internal beams will cause a change in the gain of the 
system due to some or all of the possible nonlinear effects. In the case of returning the 
transmitted beam, it appears that the gain of one of the competing patterns is 
sufficiently reduced to eliminate it from the possible patterns, 
d. The Free Standing Auxiliary Beam Pattern 

This section presents an occasionally observed result due to a small crystal 
motion. Here the crystal was translated at right angles to the beam path while the 
internal beam pattern was carefully observed for any changes. 

Note the beam configuration seen in Figure A. 13. When the crystal was 
shifted slightly to the left, the auxiliary beam pattern moved with the crystal as shown 
in Figure A. 14 instead of changing to adapt to the new crystal orientation. Note, 
however, that some slight distortion of the pattern did occur. In this case the phase- 
conjugate reflectivity was originally 32% with small amplitude oscillations on the order 
of 5°/). There was no evidence of frequency shifting of the phase-conjugate beam. 
Upon translation of the crystal, the phase-conjugate reflectivity diminished to 23% 
with oscillations of the same amplitude as before the translation. 

After translation, any rotational motion of the crystal destroys the auxiliary 
beam and a new one appears that is similar to the original pattern (i.e. Figure A. 13). 
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Figure A. 13 BaTiO^ prior to translation. 



Also, any translation in the opposite direction, causing the input beam to move over 
the auxihary beam, will destroy the auxiliary beam and a new pattern will form similar 
to the original. If the input beam is momentarily interrupted and then restored, the 
configuration appearing prior to the interruption returned. 

It is important to note that we did not always observe motion-invariance for 
patterns similar to that of Figure A. 13. The original pattern as shown in Figure A. 13 
was a commonly occurring one that could be arrived at with any of the beam profiles. 
The angle of the incident beam was not critical and neither was the choice of entrance 
face provided it was parallel to the c-axis. The pattern is usually stable and shows 
phase-conjugate reflectivity as high as 50% at some angles. 
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Figure A. 14 Free standing beams created by translating the 
cr>'stal 0.5mm from the configuration of Figure A. 13. 
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APPENDIX B 
PSIM LISTING 



PSIM_3 

A SIMULATION OF INTERNAL BEAM RESPONSE IN BARIUM TITANATE 



! IBEAM = SIZE OF INPUT BEAM 

! IBEAMF = ACTUAL NUMBER OF POINTS IN WAVEFRONT AFTER 
! DESIGNATION OF ANGLE PHI 

ISIZE = SIZE OF CRYSTAL 
ICRYS 0 = CRYSTAL GRID 

KEEP() = ARRAY TO KEEP TRACK OF CURRENTLY OCCUPIED SPACES IN THE 
GRID, 2=COL, l=ROW 
IPHI = INPUT ANGLE OF INCIDENCE 
PHI = ANGLE OF INCIDENT BEAM (REAL) 

PHIC = ANGLE OF INCIDENT BEAM AFTER ENTERING CRYSTAL 
IX = NUMBER OF COLUMN POINTS ACROSS WAVEFRONT 
lY = NUMBER OF ROW POINTS ACROSS WAVEFRONT 
ADD() = AMOUNT TO ADD TO Y FOR EACH X MOVED 
I ROW = ROW NUMBER 
! ICOL = COLUMN NUMBER 

! ITOP = HIGHEST INDEX OF REFRACTION 
I IBOT = LOWEST INDEX OF REFRACTION 

! GRAD = GRADIENT OF INDEX OF REFRACTION ACROSS THE BEAM 
! CRYS = REAL VALUE OF INDEX OF REFRACTION (BEFORE FIXING) 

! INDEX = INDEX OF REFRACTION OF LOCATION 

! REALX = PLACE FOR KEEPING TRACK OF AMOUNT TO ADD TO COLUMN 

I {OR ROW FOR REALY) WHEN MOVING DOWN THE BEAM 

! THIS ENSURES THAT ALL FRACTIONS OF A GRID SPACE 

! ARE NOT LOST 

! REALY = SEE REALX 

! REALAD () = SAME AS REALX ABOVE BUT ONE SPACE FOR EACH POINT. 

! USED DURING ACTUAL WAVEFRONT MOVEMENT 

! IPERMCO = STARTING COLUMN FOR WAVEFRONT POINT 

! IMARKB =MARKER TO KNOW WHEN REFLECTION OFF OF BOTTOM HAS OCCURRED 
! IMARKS = MARKER TO KNOW WHEN REFLECTION OFF OF SIDE HAS OCCURRED 
! ICHECK = INDEX OF TIME PERIODS 

! RSTOP = NUMBER OF TIME PERIODS IN SIMULATION 

! IDIR = NUMBER OF TIME PERIODS BETWEEN COMPUTATIONS 
! OF DIRECTION OF WAVEFRONT 

! MDIR = COUNTER FOR CHECKING AGAINST IDIR 
! XX = HORIZONTAL DISTANCE BETWEEN END POINTS OF WAVEFRONT 
! YY = VERTICAL DISTANCE BETWEEN END POINTS OF WAVEFRONT 
: ICHNG = VARIABLE USED TO CHANGE DIRECTION AFTER REFLECTION 

! ITRACO = MATRIX USED TO MARK THE INCIDENT BEAM PATH SO THAT 

! THE INDICES SET UP BY THE INCIDENT BEAM DOMINATE 

! MEAN = MEAN VALUE OF THE INDEX OF REFRACTION 
! REF = ACTUAL INDEX OF REFRACTION OF THE PHOTOREFRACTOR BEING 
! SIMULATED 

! IPLACE = POINT THAT BEAM STARTS IN CRYSTAL 

! IGRAPH = MARKER FOR POINTS TO BE PLOTTED (SET FOR 100 POINTS) 

! IGRMAX = NUMBER OF ITERATIONS BETWEEN PLOTTING WAVEFRONTS 

! CAXIS = FRACTION OF TOTAL INDEX CHANGE DUE TO ANGLE BEAM 
! MAKES WITH THEE C-AXIS IN THE CRYSTAL 

! CINDEX = REAL VALUE OF INDEX OF REFRACTION USED TO COMPUTE 
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! THE CHANGE DUE TO THE ANGLE TO C-AXIS 

! IOUT() = MATRIX FOR HOLDING VALUES OF PLOTTED POINTS 

1 ICOUT = COUNTER FOR PLACE OF POINT IN OUTPUT FILE 

! ISKIP = MARKER TO CHECK IF WAVEFRONT IS UNDERGOING REFLECTION 
! REFLCHKB =MARKER TO KNOW WHEN REFLECTION OF WAVEFRONT STARTS/ENDS 
! REFLCHKS = SAME BUT FOR SIDE 

! OPl = NUMBER OF POINTS IN WAVEFRONT TO LOOK AT FOR RECALCULATION 
! THIS EQUALS IBEAMF/PIECE 

! ISETMARK = MARKER TO TELL WHEN TO STOP LOOKING FOR REFLECTION 

! PIECE = NUMBER OF PIECES OF WAVEFRONT 

! ADDY() = FRACTION OF ONE BOX THAT EACH POINT STARTS OUT AHEAD OF 



! 

INTEGER Ibeam, Ibeamf , Isize , Icrys ( 500 , 500) 

INTEGER Keep(2.50) ,Iphi,Ix,Iy,Irow,Icol,Itop,Ibot,Reflchks(50) 
INTEGER Reflchkb(50) , Ikfac , This , That 

INTEGER Imdex, Imaexp( 50) ,Ipermc(50) ,Imarkb(50) ,Itime(50) ,Iskipb 
INTEGER Imarks (50) , lair ,Mdir , Ichng,Icout, Iskips , I temp, Opl , Piece 
INTEGER I trac ( 500 , 500) ,Mean , Iplace , Igraph . lout ( 1 , 50,101) , Isetmark 
REAL Phi,Phic ,Add(50) ,Grad,Crys,Realx,Realad(50) ,Rstop,Grmax 
REAL Xx,Yy ,Ref ,Caxis(50) , Cindex,Rchecfc,Kfac ,Realy ,Addy(50) 

I 



NOTE THAT IF BEAM SIZE EXCEEDS 50 YOU MUST REDIMENSION 



! 

ON ERROR GOSUB 5430 ! WRITES FILE TO DISK IF ERROR 

PRINTER IS 1 ! SET TERMINAL OUTPUT 

GRAPHICS OFF 

OUTPUT KBD; K; ! CLEARS SCREEN 

ON KEY 2 LABEL SAVE &STOP GOSUB 5430 
PRINT TO STOP PROGRAM AND SAVE DATA PRESS f2 

! 

I 

! 

! INPUT FACTORS 
! 

I 

Input name of data file to be used ,Name$ 

PRINT DATA FILE: ,Name$ 

INPUT SIZE OF CRYSTAL TO BE USED , Isize 
PRINT CRYSTAL SIZE IS , Isize, ON A SIDE 
INPUT BEAM SIZE IS , Ibeam 
PRINT WITH A BEAM OF SIZE , Ibeam 

INPUT PLACE ON CRYSTAL FOR BEAM TO ENTER IS , Iplace 

PRINT ENTERING AT , Iplace 

INPUT HIGH VALUE OF INDEX IS , I top 

PRINT LOW VALUE OF INDEX IS , Ibot 

PRINT INDICES RANGE FROM ,Itop, TO ,Ibot 

INPUT INITIAL ANGLE IN DEGREES IS (89 DEG MAX) , Iphi 

T1=TIMEDATE 

! 



I 



! INITIALIZE VARIABLES 

I 



Rcheck = 0 
Mdir=0 
Ref =2. 4 
lcout=0 
lsetmark=0 

Piece=5 ! REM TO CHANGE KK IF YOU CHANGE PIECE 
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1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1380 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 



PRINT NUMBER OF PIECES IN WAVEFRONT , Piece 
; 
i 



! PUT FACTORS IN OUTPUT FILE 

; 

iout(0,0,0)=lsize 
lout(0, 1 ,0)=lbeam 
lout (1,0, 0 )=Iplace 
lout(l,l,0)=ltop 
Iout(0,0,101)=Ibot 
• Iout(0,l,101)=Iphi 
Iout(l,0,101)=Piece 
lgraph=0 

Rstop=Itop*2 !SET TIME TOTAL PERIODS 

Rstop=Rstop*Isize !USE TWO LINES TO KEEP FROM INTEGER PROBLEMS 
PRINT TOTAL TIME = , Rstop 

Grmax=Rstop/100 ! SET FOR 100 WAVEFRONTS PLOTTED 
GOSUB 5170 !SET UP GRAPHICS DISPLAY 



! CONVERT PHI FROM DEGREES TO RADIANS AND FIND ANGLE INSIDE CRYSTAL 
Phi=Iphi 

Phi=Phi*(2*PI/360) 

Phic=(SIN(Phi)/Ref) 

Phic=ASN(Phic) 



X=Ibeam*COS(Phic) 

Y=Ibeam*SIN(Phic) 

Iy=Y 

Ix=X 

Ibeamf=Ix 

Opl=Ibeamf/Piece IISIZE OF PIECE TO LOOK AT WHEN CALCULATING DIRECTION 

Add(l)=Y/X 

Addl=Add(l) 

IF Iplace>(Isize-Ibeamf) THEN ! MAKE SURE IPLACE IS NOT TOO LARGE 
Iplace =Isize-Ibeamf 
END IF 



! TIME BETWEEN NONLINEAR CALCULATIONS 
Idir=5*Itop 



INITIALIZE ARRAY 



FOR 1=1 TO Ibeam 




NEXT I 



! CHANGE SIZE OF BEAM DEPENDENT ON THE INPUT ANGLE 



lrow=0 

Realy=0 

Icol=Iplace 



! FILE WAVEFRONT COORD 
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1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2530 

2590 

2600 

2510 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 



FOR 1=1 TO Ibeamf 

Keep(2,I)=Icol 
Re a iy=Re a ly+Add ( 1 ) 

Add(I)=Addtl) 

IF Realy>Irow+l THEN Irow=Irow+l 
Addy(I)=Realy-Irow 
Keep(l,I)=Irow 
Icol=Icol+l 
NEXT I 
! 

! 

! 

! 

I 

i *****SET UP INDEX OF REFRACTION FOR PATH OF INCIDENT BEAM **** 



! SET UP CRYSTAL TO MEAN INDEX 
! 

Mean=ABS((Itop+Ibot)/2) 

FOR 1=1 TO Isize 

FOR J=1 TO Isize 

Icrys(I , J)=Mean 
Itrac(I,J)=0 
NEXT J 
NEXT I 
! 

! 

! 

Grad=(Itop-Ibot)/2 ! INDEX OF REFRACTION GRADIENT 
FOR 1=1 TO Ibeamf 

Irow=Keep(l,I) 

Icol=Keep(2,I) 

Crys=Grad*SIN{ { (Ibeamf/2) - (I-l) )*(PI/Ibeamf ) ) 
Crys=Crys+Mean 
Icrys{Irow, Icol)=Crys 
NEXT I 



! SET UP VALUES DOWN THE BEAM PATH 

I 

FOR 1=1 TO Ibeamf 
Icol=Keep(2,I) 

Irow=Keep(l ,l) 

Imdex=Icrys (Irow, Icol) 

Imdexp(I )=Imdex 

Itemp=Icol 

Realx=0 

FOR J=Irow TO Isize 

Icrys(J, Itemp)=Imdex 
Itrac( J, Itemp)=l 
Realx=Realx+Add( I ) 

Itemp=Icol-Realx 
IF Itemp>Isize THEN Itemp=Isize 
IF ltemp<0 THEN ItemD=0 
NEXT J 
NEXT I 

I 

i 

1 

! 

i SET THE STARTING COLUMNS FOR THE WAVEFRONT AND 
! INITIALIZE FLAGS 



FOR 1=1 TO IBEAMF 

Ipe rmc ( I ) =Keep (2,1) 
ImarkbU)=0 
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2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 



Imarks(I)=0 
NEXT I 
! 

! 

I 

i 

I MOVEMENT 

f 

i 

! MAKE SURE BEAM IS STILL IN CRYSTAL 

f 

FOR 1=1 TO Ibeamf 

IF Imarkb(I)=l THEN 

IF Keep (1,1) <2 THEN GOTO 5140 
END IF 

IF Imarks(I)=l THEN 

IF Keep(2,I)>Isize-10 THEN GOTO 5140 
END IF 
NEXT I 
! 

i 

I 

raiLE Rcheck<Rstop 

i 

FOR 1=1 TO Ibeamf 
Icol=Keep(2,I) 

Irow=Keep(l,l) 

Imdex=Icrys(Irow,Icol) 

i 

! MARK REFLECTION IF NECESSARY 
! 

IF IcoK=l THEN 

IF lmarks(l)=0 THEN 
Imarks(l)=l 
Realad(I)=l 
END IF 
END IF 

IF Irow =Isize THEN 
Imarkb(I)=l 
END IF 



i MOVE IF IT IS TIME 
! 

IF Itime(I)<Imdex THEN 

Itime(I)=Itime(I)+l 

ELSE 

Itime(I)=0 

Realadd )=Realad( I )+Add(I ) 

IF lmarks(l)=0 THEN 

Icol=Ipermc(I)-Realad(I) 
IF lcol<0 THEN lcol=0 
ELSE 

Icol=Realad(I ) 

END IF 

Keep(2,I0=Icol 
IF Imarkab(I)=0 THEN 
Keep(l,I)=Irow+l 
ELSE 

Keep(l ,I)=Irow-l 
END IF 
END IF 
NEXT I 
! 

I 

i 
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3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3830 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 



PUT POINTS IN FILE AND DISPLAY GRAPHICS 
IF IT IS TIME TO DO SO 



Igraph=Igraph+l 
IF Igrapn>Grmax THEN 
Icout=Icout+l 
FOR L=1 TO Ibeamf 
lgraph=0 

Iout(0,L,Icout)=Keep(2,L) !X COORD OF POINT 

Iout(l,L,Icout)=Isize-Keep(l,L) !Y COORD OF POINT 
PENUP 

PLOT Iout(0,L,Icout) ,Iout(l,L,Icout) 

NEXT L 
END IF 



Rcheck=Rcheck+l 

! CHECK IF TIME TO RECOMPUTE DIRECTION ****************** 

MDIR=Mdir+l 
IF Mdir>=Idir THEN 
Mdir=0 

PRINT TABXY(1,11), CURRENT TIME PERIOD = , Rcheck 



FOR K=0 TO Piece- 1 
! 

I SET FOR 5 PIECES 

i 

IF K=0 THEN Kk=0 
IF K=1 THEN Kk=4 
IF K=2 THEN Kk=l 
IF K=3 THEN Kk=3 
IF K=4 THEN Kk=2 
!IF K=5 THEN Kk=7 
!IF K=6 THEN Kk=3 
.'IF K=7 THEN Kk=6 
!IF K=8 THEN Kk=4 
!IF K=9 THEN Kk=5 
This=Kk*Opl+l 
That=(Kk+l)*Opl 

IF Kk=Piece-l THEN That=Ibeamf 



i FIRST CHECK TO MAKE SURE THE FRONT IS NOT PARTIALLY REFLECTED 
! 

I 

lskipb=0 

lsei:mark=0 

IF Reflchkb(K)=0 THEN 
FOR M=This TO That 
IF Imarkb(M)=l THEN 
Isetmark=l 
FOR N=This TO That 

IF Imarkb(N)=0 THEN Iskipb=l 
NEXT N 
END IF 
NEXT M 
END IF 



90 



4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 



IF Isetmark=l AND lskipb=0 THEN Reflchkb(K)=l 

lskips=0 

lsetmark=0 

IF Reflchks(K)=0 THEN 
FOR M=This TO That 
IF Imarks(M)=l THEN 
FOR N=This TO That 
Isetmark=l 

IF Imarks(N)=0 THEN Iskips=l 
NEXT N 
END IF 
NEXT M 

I 

i 



IF Isetmark=l AND lskips=0 THEN Reflchks(K)=l !SKIP THIS IF NOT NECESSARY 

I 

END IF 



i 

IF lskips=0 AND lskipb=0 THEN !IF WAVEFRONT IS NOT PARTIALLY REFLECTED 



! SLOPE OF PIECE OF WAVEFRONT IS YY/XX 
I 

I 

Yy=(Keep(l ,This)+Addy(This) )- (Keep(l ,That)+Addy(That) ) 

Yl=I time ( This )/Icrys( Keep (1, This) , Keep (2, This)) 

Y2=I time ( That )/ Ic rys ( Keep(l , That) , Keep (2, That ) ) 

Yy=Yy+Yl-Y2 
IF Imarks(This)=0 THEN 

Xx=(Ipermc(That)-Realad(That) )-(Ipermc(This)-Realad(This) ) 

ELSE 

Xx=Realad(That) -Realad(This ) 

END IF 

IF Xx>l THEN 
FOR G=This TO That 

Add(G)=ABS(Yy/Xx) ! FIGURE NEW DIRECTION 
Caxis (G)=Xx/Opl 

IF Kk=Piece-l THEN Caxis(G)=Xx/ (Ibeamf-(K*Opl) ) 

NEXT G 

PRINT Add(G-l) 

END IF 
I 



FOR I=This TO That 
Icol=Keep(2,I) 
Irow=Keep(l,l) 
Imdex=Imdexp ( I ) 

! 

Itemp=Icol 

Realx=0 



! CHECK FOR REFLECTION 
! 

IF Imarkb(I)=0 THEN 

FOR J=Irow TO I size 

Icrys( J,Itemp)=Imdex 
Realx=Realx+Add( I ) 

I temp=Icol-Realx 
IF ltemp<0 THEN ltemp=0 
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4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 

5370 

3330 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

5510 



IF Itemp>Isize THEN Itemp=Isize 
Itrac( J, Itemp)=l 
NEXT J 
! 



ELSE ! INCIDENT BEAM DOMINATES 
FOR J=0 TO I row 

Ichng=Irow-J 

IF Itrac(Ichang, ltemp)=0 THEN 
Icrys{lchng,Itemp)=Imdex 
ELSE 

Icrys(Ichng, Itemp)=Icrys(Ichng, Itemp) 
END IF 
! 

Realx=Realx-Add(l ) 

Itemp=Icol-Realx 
IF ltemp<0 THEN ltemp=0 
IF Itemp>Isize THEN Itemp=Isize 
NEXT J 
END IF 
NEXT I 
END IF 
NEXT K 
END IF 

i 

1 CHECK FOR END OF RUN ******************** 

END WHILE 
T2=TIMEDATE 

PRINT TOTAL RUN TIME = ,T2-T1 
GOSUB 5360 ! WRITE FILE TO DISK 

STOP 

! SUBROUTINE TO INITIALIZE GRAPHICS 

GRAPHICS ON 

GCLEAR 

SHOW 0,lsize,0,lsize 
PLOT 0,0, -2 
FOR 0=0 TO I size 
PLOT 0,0, -1 
NEXT 6 

FOR Z=0 TO Isize 
PLOT Isize, Z,-l 
NEXT Z 

FOR Q=Isize TO 0 STEP -1 
PLOT Q, Isize, -1 
NEXT Q 

FOR Z=Isize TO 0 STEP -1 

PLOT 0,Z,-1 

NEXT Z 

PLOT 0,0, -2 

RETURN 

! SUBROUTINE TO OUTPUT FILE TO DISK 
MASS STORAGE IS :, 700,0 
CREATE BOAT Name$ , 10404,2 
ASSIGN ?ath_l TO Name 3 
OUTPUT ?ach_l;Iout(^) 

ASSIGN ?aT:h_l TO ^ 

RETURN 

! SUBROUTINE TO OUTPUT FILE IF ERROR OCCURS 
MASS STORAGE IS :, 700,0 
CREATE BDAT Name$ , 10404,2 
ASSIGN Path_l TO Name$ 

OUTPUT Path_l;Iout(*) 

ASSIGN Path_l TO * 

PRINT ERROR! PROGRAM TERMINATED 
T2=TIMEDATE 

PRINT TOTAL RUN TIME = ,T2-T1 
STOP 
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5520 RETURN 
5530 end 
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